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Abstract 

The COMPASS collaboration has collected the currently largest data set on diffractively 
produced 7r“7r“7r+ final states using a negative pion beam of 190 GeV/c momentum impinging 
on a stationary proton target. This data set allows for a systematic partial-wave analysis 
in 100 bins of three-pion mass, 0.5 < tostt < 2.5GeV/c^, and in 11 bins of the reduced four- 
momentum transfer squared, 0.1 < t' < 1.0(GeV/c)^. This two-dimensional analysis offers 
sensitivity to genuine one-step resonance production, i.e. the production of a state followed 
by its decay, as well as to more complex dynamical effects in nonresonant Stt production. 
In this paper, we present detailed studies on selected Stt partial waves with = 0 
1++, 2 2++, and 4++. In these waves, we observe the well-known ground-state mesons 

as well as a new narrow axial-vector meson ai(1420) decaying into /o(980)7r. In addition, 
we present the results of a novel method to extract the amplitude of the 7r“7r+ subsystem 
with = 0+ O'*'’*' in various partial waves from the 7r“7r“7r+ data. Evidence is found 

for correlation of the /o(980) and /o(1500) appearing as intermediate 7r“7r+ isobars in the 
decay of the known 7r(1800) and 7r2(1880). 
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1 Introduction 


In this paper, we report on the results of a partial-wave analysis of the 7r“7r“7r'’' system produced 
by a 190GeV/c 7r“ beam impinging on a liquid-hydrogen target. The reaction of interest is 
diffractive dissociation of a 7r“ into a 7r“7r“7r+ system, 

TT~ + P^ 7r“7r“7r+ -h Precoil, (1) 


with Precoil denoting the recoiling target proton. The data for this analysis were recorded with 
the COMPASS experiment at the CERN SPS in 2008. 


Despite many decades of research in hadron spectroscopy, the excitation spectrum of light mesons, 
which are made of u, d, and s quarks, is still only partially known. In the framework of the 
simple constituent-quark model using SU(3)gg^.^^Qj. (8> SU(2)j,pjj^ (g) SU(3)pjj[Qj, symmetry, a number 
of frequently observed states are commonly interpreted in terms of orbital and radial excitations 
of quark-antiquark ground-state mesons, i.e. they are assigned to the multiplets resulting from 
the symmetry. Some of these assignments are still disputed, as e.g. the isovector mesons p(1450), 
p(1700), 7r(1300), and 7r(1800) [J, as well as the whole sector of scalar mesons j2]. In addition, 
a number of extra states have been found, which cannot be accommodated by the constituent- 
quark model. These extra states appear in mass ranges where quark-model states have already 
been identified, e.g. the 7r2(1880) which is close to the 7r2(1670) ground state. Other observed 
states seem to have peculiar decay modes or decay widths that do not fit well into the general 
pattern. Searching for new states beyond the constituent-quark model, attempts have been made 
to establish the existence of gluonic deg rees of freedom. The fingerprints are expected to be 
so-called exotic spin quantum number! I®'! or decay branching ratios, which could identify them as 
hybrids Ell], glueballs mi, or tetra-quark systems jl]. Potential candidates are e.g. 7ri(1600), 
7r(1800), 7r2(1880) or /o(1500), /o(1710) or /o(980), ao(980), /i(1420), respectively. 


The COMPASS collaboration has already studied properties of isovector 37r resonances 018] in 
the mass range between 1.1 and 2.1 GeV/c^ using a lead target. In this paper, isovector mesons 
decaying into three charged pions are studied using a hydrogen target with the emphasis on i) 
production kinematics, ii) separation of nonresonant processes. Hi) search for new and excited 
mesons, and iv) on properties of the vrvr 5-wave amplitude. This paper is the hrst in a planned 
series of publications to present precision studies revisiting all quantum numbers accessible in 
reaction Q up to total spin J = 6. The analysis is limited to states belonging to the family of vrj 
and oj. In addition, the large data set allows to apply a novel method for investigating isoscalar 
states, which occur as TT~Tr~^ subsystems in the decays of isovector mesons. 


The Particle Data Group (PDG) |9| lists a total of eleven well-established isovector states with 
masses below 2.1 GeV/c^ (see Table Q, where only the oq states do not decay into Sir due to 
parity conservation. The widths of the oq, 0 , 2 , and 04 ground states have values of about 10% of 
their mass values, while the ai(1260) is much broader. Pionic excitations are typically broader 
with values of their width being about 15 to 20 % of their mass values. In addition, the table 
contains a number of less well-established states. Even for some established states, properties 
such as mass and width are poorly determined, e.g. for the ai(1260) as the lightest oi state, 
the reported widths vary between 250 and 600 MeV/. Another example is the inconsistency 
in the mass measurements of 7r(1800), where experimental results cluster around two different 
mean values. This has lead to speculations on the existence of two states, one being an ordinary 
meson and the other one a hybrid. Extensive discussions of the light-meson sector are found in 
Refs. HI ITO] . 

The partial-wave analysis of the Stt system has a long history [T]. The technique of partial-wave 


|a| jPC numbers that are forbidden for qq in the nonrelativistic limit. 
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Table 1: Resonance parameters of aj and ttj mesons in the mass region below 2.1 GeV/c^ as given in 
PDG [S]. Note that due to parity conservation the Oq states cannot decay into 7r“7r“7r+. 


Particle 

jPC 

Mass [MeV/c2] 

Width [MeV/c2] 



Established states 

ao(980) 

0 ++ 

980 ± 20 

50 to 100 

ai(1260) 

1 ++ 

1230 ± 40 

250 to 600 

02(1320) 

2 ++ 

1318.3t;);^ 

107 ±5 

ao(1450) 

0 ++ 

1474 ± 19 

265 ± 13 

04(2040) 

4++ 

1996lg° 

255128 

7r(1300) 

0 -+ 

1300 ± 100 

200 to 600 

7ri(1400) 

1-+ 

1354 ± 25 

330 ± 35 

7ri(1600) 

1-+ 

16621® 

241 ± 40 

7r2(1670) 

2 "+ 

1672.2 ±3.0 

260 ±9 

7r(1800) 

0 -+ 

1812 ± 12 

208 ± 12 

712(1880) 

2 "+ 

1895 ±16 

235 ± 34 


States omitted from summary table 


ai(1640) 

1++ 

1647 ± 22 

254 ± 27 

02(1700) 

2++ 

1732 ±16 

194 ± 40 

7r2(2100) 

2"+ 

2090 ± 29 

625 ± 50 



Further states 


03(1875) 

3++ 

1874 ± 43 ± 96 

385 ±121 ±114 

oi(1930) 

1++ 

19301?!) 

155 ± 45 

02(1950) 

2++ 

19501?!) 

180;?!) 

02(1990) 

2++ 

2050 ± 10 ± 40 

2003 ± 10 ± 19 

190 ± 22 ± 100 

249 ± 23 ± 32 

oo(2020) 

0++ 

2025 ± 30 

330 ± 75 

02(2030) 

2++ 

2030 ± 20 

205 ± 30 

03(2030) 

3++ 

2031 ± 12 

150 ± 18 

oi(2095) 

1++ 

2096 ± 17 ± 121 

451 ±41 ±81 

712(2005) 

2"+ 

1974 ± 14 ± 83 

2005 ±15 

341 ± 61 ± 139 

200 ± 40 

7ri(2015) 

1-+ 

2014 ± 20 ± 16 

2001 ± 30 ± 92 

230 ± 32 ± 73 

333 ± 52 ± 49 

7r(2070) 

0-+ 

2070 ± 35 

310 ;?°° 

X(1775) 

7- + 

1763 ± 20 

1787 ± 18 

192 ± 60 

118 ±60 

X(2000) 

?■ + 

1964 ± 35 
~ 2100 

2214 ± 15 

2080 ± 40 

225 ± 50 
~ 500 

355 ± 21 

340 ± 80 












5 


analysis (PWA) of Svr systems was established by the work of Ascoli et al. jlll [12] in 1968. The 
CERN-Munich collaboration (ACCMOR) [131 O CSl CS] further developed this method and 
measured signihcant contributions from partial waves up to J = 2, without including spin-exotic 
waves. The largest data set used so far, which is the basis of several publications on the 3tt hnal 
state, was obtained and analyzed by the BNL E852 collaboration nziiiiKini. They have studied 
reaction Q at beam momenta of 18GeV/c and observed signihcant waves with = 0 
1^"*“, 2~^^, and 2 quantum numbers. In addition, they have detected a 1 spin-exotic wave 
in the /?(770) vr decay channel with signihcant huctuation in intensity depending on the number 
of partial waves used, i.e. with a considerable model dependence. Also the VES experiment 
has large data sets, the analysis of which was published mostly in conference proceedings, see 
e.g. Refs. [2011211 [221123]. 



TT' 


TT' 


TT 


+ 


Precoil 


Figure 1: Diffractive dissociation of a beam pion on a target proton into the three-pion final state. The 
figure shows the excitation of an intermediate resonance X~ via Pomeron exchange and its subsequent 
decay into Stt. 


As illustrated in Fig. for reaction 0 at 190 GeV beam energy, the strong interaction can be 
described by the exchange of a quasi-particle called Pomeron, P, which is a fiavorless glueball-like 
object that accounts for diffractive dissociation and most of the two-body elastic scattering |24j . 
The Regge trajectory ap{t) of the Pomeron determines the elastic scattering amplitude 

A(s,t) oc (2) 


Here, s is the squared center-of-mass energy, t the squared four-momentum transferred between 
beam particle and target nucleon, and 


Oip{t) — 1 -|- 6 p CTp t, (3) 

where 0.081 ^ ep < 0.112 and Up ~ 0.25 (GeV/c)“^ [24] . The Pomeron is an even-signature 
Regge trajectory with ..., and its first Regge pole corresponds to a 

flavorless hadron with = 2^~^ and a mass of about 1.9GeV/c^. The parameter a'p modifies 
the dependence of the differential cross section on the four-momentum transfer. Equation Q 
implies a dependence of the cross section on t as 


The slope parameter b is given by 


dcr 

dt 


oc e 


—b t 


(4) 


b 


bo + Aa'p In 



(5) 


where bo is a generic slope parameter and the unknown scale parameter sq is usually taken to be 
1 GeV^. The reduced four-momentum transfer squared is 


t' = |t| - |t|min > 0, where |t|min 


f ml^-m 

y 2|pbeam| 


2 

TT 


2 


( 6 ) 
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Table 2: List of allowed quantum numbers for X assuming that it is produced in the interaction of 
a = 0 beam pion and a 2'^'^ Pomeron as an example, with relative orbital angular momentum ^ 
between the two. 


1 

JPC 

of X 



0 

2 -+ 





1 

1 ++, 

2 ++, 

3++ 



2 

+ 

1 

o 

1-+, 

2 -+, 

CO 

1 

+ 

+ 

1 

CO 

1 ++, 

2 ++, 

3++, 

4++, 

5++ 

4 

2 -+, 

1 

CO 

+ 

1 

5-+, 

6“+ 


is the minimum momentum transfer needed to excite the beam particle to a mass msjr, which 
is the invariant mass of the Stt hnal state. The beam momentum pbeam is measured in the 
laboratory frame. For the Svr mass range of 0.5 to 2.5GeV/c^ considered in this analysis, typical 
values of |t|min are well below 10“^ (GeV/c)^. Different production mechanisms, i.e. different 
exchange particles, can lead to different slopes b. The existence of concurrent exchange processes 
thus results in a more complex form of the t' dependence with coherently and/or incoherently 
overlapping exponentials. The t' range for this analysis is 0.1 to 1.0 (GeV/c)^. 

Studies of diffractive dissociation of pions, see e.g. Refs. [g dg [221 ES], reveal the existence of 
nonresonant background processes such as the Deck effect |25j . These processes exhibit strongly 
mass-dependent production amplitudes that occur in the same partial waves as the resonances 
under study. In particular, the analyses presented in Refs. |16l fig showed the importance of the 
kinematic variable t' in a partial-wave analysis and illustrated the power of accounting for the 
difference in the t' dependence of the reaction mechanisms and also of the different resonances. 
In this work, we take advantage of the large size of our data sample and develop this approach 
further in order to better disentangle resonant and nonresonant components. 


In the case of Pomeron exchange, the partial waves induced by a pion beam can be assessed as 
follows: the 7r“ is an isovector pseudoscalar with negative G parity and the Pomeron is assumed 
to be an isoscalar C = -|-1 object, so that the partial waves all have = 1“. Possible 
quantum numbers W of partial waves are listed in Table for the lowest values of the relative 
orbital angular momentum i between the beam particle and a Pomeron as an example. 

As we will demonstrate in this paper, almost all partial waves listed in Table are indeed observed 
in our data. Higher-spin waves with J > 5 contribute significantly only at masses above 2 GeV /(?. 
The table includes spin-exotic partial waves such as = 1 3 and 5 The present 

paper focuses on non-exotic spin quantum numbers with the emphasis on known states. They are 
extracted from the data by partial-wave methods that contain an a priori unknown dependence 
on t', which is extracted from the data. 


The work related to this topic is subdivided into two publications, owing to the large amount of 
material and various, in parts novel analysis techniques used. This paper contains details on the 


experiment in Section 2.1 and a description of the basic event selection criteria in Sections 2.2 


and 2.3 where we also present the general features of our data set and the overall kinematic 
distributions for both msTr and t'. Section contains a detailed description of our analysis 
method and the PWA model used. For clarity, we include a rather extensive mathematical 


Although the C parity is not defined for a charged system, it is customary to quote the quantum numbers 
of the corresponding neutral partner state in the isospin multiplet. The C parity can be generalized to the G parity 
G = a multiplicative quantum number, which is defined for the non-strange states of a meson multiplet. 
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SciFi and 
beam counter 


Non-interacting 

beam 



hodoscopes barrels veto veto 


Figure 2: (color online) Simplified scheme of the diffractive trigger. The main components are the beam 
trigger, which selects beam particles, and the RPD, which triggers on slow charged particles leaving the 
target. The veto system (red) rejects uninteresting events and consists of three parts: The veto hodoscopes, 
the sandwich, and the beam veto. 


description summarizing the work of many authors, who laid the basis for our analysis (see 
e.g. Refs. [26l [271 EH EH ISQl [SH [32]) • In this scheme, the analysis follows a two-step procedure 
described in Ref. |32j . In the first step, a PWA is performed in bins of and t'. The results 
of this so-called mass-independent fit are presented and discussed in Sections and In these 
and the following sections, the focus lies on 37r resonances with masses below 2.1 GeV/c^. The 
discussion on t' dependences includes the kinematic distributions and J^'^-resolved t' spectra. 
In Section]^ we present a novel approach that allows us to investigate the amplitude of Tr^Tr"*" 
subsystems in the decay process. In particular, we address the topic of the scalar sector containing 
/o mesons and its complicated relation to tt-tt S-wave scattering. The relation of /o(980) and 
/o(1500) mesons to vrTr scattering will be demonstrated. In this paper, all error bars shown in 


the figures represent statistical uncertainties only. Systematic effects are discussed in Section 4.6 


and Appendix [B| In Section we conclude by summarizing the hndings based on qualitative 
arguments. The appendices contain details about more technical issues. 


The analysis methods and results presented in this paper will serve as a basis for further 
publications that will be dedicated to individual partial waves. In the second step of the analysis, 
physics parameters will be extracted from the data presented in this paper by performing a ht 
that models the resonance amplitudes and the amplitudes of nonresonant processes. This involves 
simultaneous htting to many partial-wave amplitudes in all bins of t'. Such a mass-dependent 
ht, which will allow us to extract the t' dependences of various components, i.e. resonant and 
nonresonant contributions for individual partial waves as well as resonance parameters for the 
mesonic states observed with different , will be described in a forthcoming paper |33j . 


2 Experimental Setup and Event Selection 
2.1 COMPASS Setup 

The COMPASS spectrometer, which is described in general in Ref. |34| . is situated at the CERN 
SPS. The setup used for the measurement presented here is explained in more detail in Ref. |35) . 
COMPASS uses secondary hadron and tertiary muon beams that are produced by the 400GeV/c 
SPS proton beam impinging on a 50 cm long beryllium target. The measurement described in 
this paper is based on data recorded during the 2008 COMPASS run. The beam was tuned to 
deliver negatively charged hadrons of 190 GeV/c momentum passing through a pair of beam 
Cherenkov detectors (CEDARs) for beam particle identihcation. The beam impinged on a 40 cm 
long liquid-hydrogen target with an intensity of 5 x 10^ particles per SPS spill (10 s extraction 
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with a repetition time of 45 s). At the target, the hadronic component of the beam consisted of 
96.8% 7r“, 2.4% K~, and 0.8% p. In addition, the beam contained about 1 % //“ and an even 
smaller amount of electrons. 

The target was surrounded by a Recoil-Proton Detector (RPD) consisting of two concentric, 
inner and outer, barrels of scintillators with 12 and 24 azimuthal segments, respectively. Recoil 
protons emerging from diffraction-like reactions must carry momenta of at least 270 MeV/c in 
order to traverse the target containment and to be detected in the two RPD rings. This leads to 
a minimum detectable squared four-momentum transfer t' of about 0.07 (GeV/c)^. 

Incoming beam particles and outgoing reaction products that emerge in the forward region were 
detected by a set of silicon micro-strip detector stations, each consisting of two double-sided 
detector modules that were arranged to view four projections. Particles emerging in the forward 
direction were momentum-analyzed by the two-stage magnetic spectrometer with a wide angular 
acceptance of ±180mrad. Both spectrometer stages are each composed of a bending magnet, 
charged-particle tracking, electromagnetic and hadronic calorimetry, and muon identification. 
Particles in the momentum range between 2.5 and 50GeV/c and passing through the Ring- 
Imaging Cherenkov (RICH) detector in the first stage can be identified as pion, kaon, or proton. 
The experiment offers large acceptance and high reconstruction efficiency over a wide range of 
three-pion mass and squared four-momentum transfer t'. 

2.2 Hardware Trigger 

A minimum-bias trigger, the so-called diffractive trigger (DTO) |35[ I36j . was used to preselect 
events with interacting beam particles and a recoiling proton emerging from the target. The 
trigger elements are shown schematically in Fig. The DTO trigger is a coincidence of three 
independent trigger signals: i) the beam trigger, ii) the recoil-proton trigger, and Hi) the veto 
signal. Incoming beam particles are selected by the beam trigger requiring a signal in one plane 
of the scintillating-fiber detector (SciFi) in coincidence with a hit in the beam counter, which 
is a scintillator disc of 32 mm diameter and 4 mm thickness. Both beam-trigger elements are 
located upstream of the target. The proton trigger selects events with protons recoiling from the 
target. It features target pointing and discrimination of protons from other charged particles by 
measuring the energy loss in each ring of the RPD. The veto signal has three sub-components. 
The veto hodoscopes reject incoming beam particles with trajectories far from the nominal one. 
Similarly, the sandwich scintillation detector that is positioned downstream close to the target, 
vetoes events with particles leaving the target area outside of the geometrical acceptance of the 
spectrometer. Lastly, the beam veto, two scintillator discs of 35 mm diameter and 5 mm thickness 
positioned between the second analyzing magnet and the second electromagnetic calorimeter, 
vetoes signals from non-interacting beam particles. Events recorded with the diffractive trigger 
can be regarded as good candidates for diffractive dissociation reactions. 


2.3 Event Selection 


The analysis is based on a data set of about 6.4 x 10® events selected by the hardware trigger (see 
Section 2.2). The event selection aims at a clean sample of exclusive Tr~ -|-p —)• 7r“7r“7r+ + Precoil 
events (see Fig. Q and consists of the following criteria (see Ref. m for more details): 


i) A vertex is required to be formed by the beam particle and three charged outgoing tracks 
with a total charge sum of —1. The vertex must be located within the fiducial volume of 
the liquid-hydrogen target. 

ii) Momentum conservation is applied by requiring exactly one recoil particle detected in the 
RPD that is back-to-back with the outgoing 7r“7r“7r'’“ system in the plane transverse to the 
beam (transverse momentum balance). This suppresses contributions from double-diffractive 




9 


xlO® 



Figure 3: Panels (a) and (b) show the reconstructed beam energy Fbeam after selection cuts (filled 
histograms). The open histogram in (a) represents the energy distribution without the RPD information. 
In the zoomed view (b), the vertical red lines indicate the accepted range. 


processes, in which also the target proton is excited. 
iii) The beam energy i^beamj which is calcnlated from the energy and momentnm of the three 
outgoing particles corrected for the target recoil, must be within a window of ±3.78 GeV 
around the nominal beam energy, which corresponds to two standard deviations (see Fig. 3b). 


A number of additional selection criteria is applied in order to reject background events originating 
from other processes. Events are disregarded if the incoming beam particle is identified by the two 
beam Cherenkov detectors (CEDARs) as a kaon. This suppresses kaon-beam induced events, like 
e.g. K~ ±p —)• ±Precoil- If Sit least one of the three forward-going particles is identified 

by the RICH detector as a kaon, proton, electron, muon, or noise, the event is rejected, thereby 
suppressing events such as e.g. vr” ±p —>■ tt~K~K^ TPrecoil- In order to reject background events 
stemming from the central-production reaction vr” ± p —)• ± 7r“7r'*' ± Precoil) in which no 

three-pion resonances are formed, the faster tt~ in the event is required to have a Feynman-x 
below 0.9 defined in the overall center-of-mass frame. The rapidity difference between the faster 
7 r~ and the remaining vr^Tr"’' pair is limited to the range from 2.7 to 4.5. Figure]^ shows the 
and m^-^+ distributions of the sample that is cut away. 


After all cuts, the data sample consists of 46 x 10® events in the analyzed kinematic region 
of three-pion mass, 0.5 < msTr < 2.5GeV/c^, and four-momentum transfer squared, 0.1 < t' 
< 1.0 (GeV/c)^. Figures 5a and 5b show for all selected events the mass spectrum of Tr^Tr^Tr"*" 


Figures 

and of the two 7 r~'K~^ combinations. The known pattern of resonances ai(1260), 02(1320), and 


7r2(1670) is seen in the 37r system as well as p{770), /o(980), /2(1270), and /93(1690) in the 7r“7r'’“ 
subsystem. From Fig. 5c the correlation of the resonances in the 7r“7r“7r+ system and in the 
7r“7r“’' subsystem is clearly visible. This correlation is the basis of our analysis model described in 


Section The t' spectrum is shown in Fig. 5d 


A Monte Carlo simulation has shown that for the reaction under study, the 37r mass resolution of the 
spectrometer varies between 5.4MeV/c^ at small (in the range from 0.5 to l.OGeV/c^) and 
15.5 MeV jc? at large 7713 ,^ (in the range from 2.0 to 2.5 GeV/c^), respectively. The t' resolution as 
obtained from the reconstructed 37r final state ranges between 7 x 10“^ and 20 x 10“® (GeV/c)^ 
depending on the majr and t' region. The resolution of the reconstructed beam energy Fibeam is 
smaller than the intrinsic energy spread of the beam and varies between 0.6 and 0.9 GeV. The 
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Figure 4: (color online) Effect of the central-production (CP) veto. Panel (a): 7r“7r“7r+ invariant mass 
spectrum without the central-production veto (yellow histogram) together with the sample that is removed 
by the central-production veto (red histogram). The inset shows the same histogram with magnified 
ordinate scale. Note that the partial-wave analysis is performed only in the mass region of 0.5 < 

< 2.5 GeV !(? indicated by the vertical red lines. Panel (b): 7r“7r“'' invariant mass distribution (two entries 
per event) of the sample that is cut away. 


position of the primary interaction vertex along the beam axis is reconstructed with a resolution 
of approximately 6 mm at small and 1.5 mm at large mz-K- The overall detection efficiency, 
which includes detector acceptance, reconstruction efficiency, and event selection, is estimated for 
isotropically distributed (phase-space) 7 r“ 7 r“ 7 r'’' events. Integrated over the analyzed kinematic 
region, it is 49% on average. More details are found in Appendix and Ref. |37| . 


3 Partial-Wave Analysis Method 

The goal of the analysis described in this paper is to extract the resonances contributing to the 
reaction 7 r“ -|-p —)• Tr^Tr^Tr"*" -fprecoil and to determine their quantum numbers from the observed 
kinematic distributions of the outgoing 7 r“ 7 r“ 7 r'’' system. This is accomplished using partial-wave 
analysis techniques. The basic assumption is that resonances dominate the Svr intermediate states 
X~ produced in the scattering process, so that the X~ production can be treated independently 
of the X~ decay (see Fig. [^. The amplitude for a certain intermediate state X~ is therefore 
factorized into two terms: i) the transition amplitude 'T describing the production of a state X~ 
with specific quantum numbers and ii) the decay amplitude iF that describes the decay of the 
X~ state into a particular 7 r“ 7 r“ 7 r^ final state. For fixed beam energy, the measured kinematic 
distribution of the final-state particles depends on the Svr invariant mass 772 . 37 ^, the four-momentum 
transfer squared t', and a set of five additional phase-space variables denoted as r, which fully 
describe the three-body decay and are defined below. 

3.1 Isobar Model 

In order to illustrate the isobar ansatz, we give in Fig. |^two examples for Dalitz plots for two 
different regions of 7723 , 7 . In the Svr mass region around 02(1320), which also includes contributions 
from ai(1260), we see a dominant contribution of the p(770) in the subsystem, while 

for values of 7723 ,, around 772(1670) several 27r resonances contribute, i.e. /5(770), /o(980), and 
/2(1270). 
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(b) [GeV/c2] 



Figure 5: (color online) Final event sample after all selection cuts: (a) invariant mass spectrum of 
7r“7r“7r''' in the range used in this analysis (see vertical lines in Fig. 4a), (b) mass distribution, 

(c) correlation of the two, (d) t' distribution with vertical lines indicating the range of t' values used in 
this analysis. The histograms in (b) and (c) have two entries per event. The labels indicate the position 
of major Stt and 27r resonances, the gray shaded areas in (a) the mass regions used to generate the Dalitz 
plots in Fig. 1^ 
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Figure 6: (color online) (a) Dalitz plot in the mass regions of the 02(1320), which also includes the 
ai(1260), (b) around the 7r2(1670). The used Stt mass regions are indicated in Fig. 5a The dominant 
p(770) TT decays of ai(1260) and 02(1320) are clearly visible. The 712(1670) region exhibits p(770) tt, 
/2(1270) tt, and /o(980) tt decay modes. 
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TT 
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Figure 7: The decay oi X , as described in the isobar model, is assumed to proceed via an intermediate 
7r“7r+ state the so-called isobar. 


Because of the strong contribution of resonances in the tt vr'*' subsystem, _;^e three-body decay 
amplitude !J'(r, ms,, 


is known as the isobar modei 1'^ 


is factorized into two two-body decay terms (see Fig. 
and the introduced intermediate vr 


vr' 


This factorization 
state ^ is called the isobar. 
In the hrst two-body decay, X~ —vr”, a relative orbital angular momentum L is involved 
in the decay. The decay amplitude ?F(r, msTr) completely describes the kinematic distribution 
of the three outgoing pions for particular quantum numbers of X~ and for a particular isobar 
channel with a given L. 


The two subsequent two-body decays are described in different right-handed coordinate systems, 
i.e. the Gottfried-Jackson and the helicity reference frame (see Fig.[^. The Gottfried-Jackson (GJ) 
frame is used to describe the angular distribution of the decay of the intermediate state X~ 
into the isobar ^ and the bachelor pion. It is constructed in the X~ rest system, in which the 
direction of the beam particle defines the zqj axis and the yoj is oriented along the normal 
to the production plane (yoj = Pbeam ^ Px'° ~ ^recoil ^ Pb^am’ where unit vectors are indicated 
by a circumflex). In this system, the momenta of the isobar and the bachelor pion are back to 
back, so that the two-body decay X~ —7r“ is described by the polar angle -dcj and the 
azimuthal angle (/)ty of the isobar, the latter being also referred to as Treiman-Yang angle. 


early detailed discussion can be found in Ref. m- 
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VGJ = Precoil X Phlm 



^Gj TT Ph. 


GJ 

earn 


X production plane 


Figure 8: (color online) Definition of the Gottfried-Jackson reference frame (GJ) in the X rest system 
and of the helicity reference frame (HF) in the rest system as they are used to analyze the angular 
distributions of the decays X~ —-I- 7r“ and 7r“ -|- 7r+, respectively. Unit vectors are indicated by 

a circumflex. 


For the decay of the isobar ^ into vr^vr"’', the helicity reference system (HF) is used to describe 
the angular distribution. This frame is constructed by boosting from the Gottfried-Jackson 
system into the ^ rest frame. The zhf axis is taken along the original direction of the isobar 
and yjjF = -^GJ x ipF- The two pions are emitted back to back, so that the —)• decay is 

described by the polar angle i?hf and the azimuthal angle (/>hf of the negative pion. 

For illustration, Fig. shows the observed, i.e. acceptance-uncorrected angular distributions in 
the two reference systems for events around the 7r2(1670) mass region. The main decay of this 
resonance is through the /2(1270) isobar, which is a = 2'^'^ state decaying into vr^Tr"’' in a 
relative D-wave in the helicity reference frame. The /2(1270) and the bachelor pion are emitted 
in a relative S or D-wave in the Gottfried-Jackson coordinate system. Note that the shown 
distribution is complicated by the fact that other decay modes of the 7r2(1670) as well as decays of 
other Stt resonances with different angular distributions interfere with the 7r2(1670) —>■ /2(1270) Tr~ 
decay. 


3.2 Parametrization of Decay Amplitudes 

In the helicity formalism |261 l?7l 131) . the amplitude for a two-body decay of a state R 
with spin J into particles 1 and 2 can be factorized into a dynamic part that 

describes the mass dependence and an angular part. The latter is related to the rotation between 
the rest frame of the parent system R, in which its spin projection M is defined, and the helicity 
frame used to define the daughter spin states, which are given by the helicities Ai^ 2 - The rotation 
is described by the Wigner ZJ-function. In addition, there are two Glebsch-Gordan coefficients 
arising in the decay Z? —)• 1 J- 2: i) for the coupling of the spins Ji ,2 of the daughter particles to 
the total intrinsic spin S and ii) for the coupling of the relative orbital angular momentum L 12 
between the daughter particles with S to J. As the orbital angular momentum L 12 in the decay 
is by definition perpendicular to the quantization axis in the helicity formalism, its z projection 
vanishes. 

The amplitude for the two-body decay of the isobar ^ with spin and helicity A into Tr^Tr"*" 
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Figure 9: (color online) Example of a Stt angular distribution observed in the mass region 1.6 < TO 3 ,r 
< 1.7GeV/c^ around the 712(1670) indicated by vertical red lines in the upper left panel. The main decay 
of this resonance is through the /2(1270) isobar, which is a = 2++ state, decaying into 7r“7r+ in a 
D-wave. The /2(1270) and the bachelor pion are in a relative S or D-wave. 


is given by 

>^i('f?HF, m^) = ('^hf, -i^HF, 0) (7) 

with being the vr^vr"’' invariant mass. The dynamic part factorizes into several components: 

ag Tj^(m^; m^) Z\g(m^; m^). (8) 

'-v-" '-V-' 

normalization dynamics 

Here, the fact was already used that pions are spinless isospin-1 particles. Therefore, the L-S 
coupling Clebsch-Gordan coefficient is unity and the orbital angular momentum in the decay 
is identical to the spin of the isobar. The coupling amplitude describes the strength of the 
decay and is usually unknown. Parametrizations for the barrier factor Fj^ and the isobar line 
shape are discussed in Section 

The amplitude for the two-body decay of X~ into the isobar ^ and the bachelor pion is 
constructed by summing over the helicity A of the intermediate isobar: 

= ^T>ma('/>ty,??gj, 0)/Ao("r37r;mg,m^). (9) 

A 

Taking into account the quantum numbers of the bachelor pion the dynamic part of the amplitude 
reads: 

fxoi'^37T',rn^,mT,) = x/2L -b 1 (L 0 A | J A) ax m^, m^). (10) 

--^ ^ ^ 

normalization L-S coupling dynamics 

Clebsch-Gordan 

This is the nonrelativistic L-S coupling scheme as introduced by Jacob and Wick in Ref. |26| . 
which is equivalent to the nonrelativistic Zemach tensors |381 [39]. Relativistic corrections as 
worked out in Ref. m are not applied. The results presented here are therefore comparable to 


4.1 
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those of previous analyses. The relativistic corrections are expected to become important for large 
breakup momenta in the X~ —?• + vr” decay and will be studied in detail in a future analysis. 


In Eq. (101, again an unknown coupling amplitude ax appears. Note that the line shape Axinis-^) 
of the X~ is unknown. It is actually the goal of the analysis to extract it from the data. This 
is achieved by setting Ax to unity so that it does not appear in the above formula and by 
performing the analysis in narrow bins of msTr, thereby neglecting the 1713 ^ dependence within 
each bin. 


Combining Eqs. Q and ([^ yields the X decay amplitude 


= r 

rn^, ^GJ, </'TY; W-Stt) = ^ D-h , '&GJ, 0) /Ao("^37r; ^^('dHF, 

A 

( 11 ) 

However, the above amplitude does not yet have the correct Bose symmetry under exchange of 
the two indistinguishable 7r“ in the hnal state. The symmetrized amplitude is 


i^ij(Tl3,'F23;m3^) 


1 

71 


[tpi,jins', msTr) + '>Pi,jins', 


( 12 ) 


where ri 3 and r 23 are the sets of phase-space variables calculated for the two possible 7r“7r''“ 
combinations of the 7rj|"7r^7r^ system. Equation (12| takes correctly into account the self¬ 
interference due to the particle-exchange symmetry. For better readability, we will use the 
simplihed notation ^ijiT^m^T^) in the text below. 


The X~ decay amplitude 'I'ij is uniquely defined by two indices: a) the set of X~ quantum 
numbers (isospin I, G parity, spin J, parity P, C parity, and the spin projection M), represented 
here by the index i = (J*^, , M), and b) by the X~ decay mode enumerated by j = (^, L). In 

this way we can describe the decay of a diffractively produced intermediate state X~ with mass 
m 37 r decaying into a 7r“7r'’“ isobar ^ and a bachelor 7r“. 


3.3 Partial-Wave Decomposition 

The intensity distribution Xim^T^^t', t) of the final-state particles is written as a truncated 
series of partial waves denoted by the indices i and j, which represent certain quantum number 
combinations as discussed in Section |3.2[ The strengths and phases, with which the different 
intermediate states X~ are produced, are described by the production amplitudes Tnim^i^^t'). 
They depend on the production kinematics and on the set i = M) of the X~ quantum 

numbers. Together with the decay amplitudes from Eq. ( |12[ ), the intensity is written as the 
coherent sum over the different intermediate X~ states represented by i and the different X~ 
decay modes enumerated by j: 


m 


e=±lr=l 




(13) 


In the above formula, two additional indices, the so-called reflectivity e and the rank index r, are 
introduced, which are both summed over incoherently. Before discussing these two indices, we 


transform Eq. (13) further. 


In the helicity formalism, the isobar-model decay amplitudes are calculable up to the unknown 

^ [see 

Eqs. 


couplings a^ and ax, which appear at each decay vertex and were introduced in Section 3.2 

and (10) and Fig. [^. Assuming that these couplings do not depend on the kinematics. 
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these unknowns can be be pulled out of the decay amplitude in Eq. (12) and absorbed into the 
production amplitudes by the following redehnitions: 






ax 


(14) 


and 


TTj (?^37r, t') = a^ ax 77^ t'). 


(15) 


Note that now the amplitudes carry not only information about the production of the state i, 
but also about its coupling to a certain decay channel j. Therefore, we will refer to the as 
transition amplitudes in the rest of the text. We introduce the index 


a={i,j). (16) 

This notation represents a certain partial wave and contains all information about the production 


as well as the decay (see Section 3.2). With these modifications, we rewrite the expression for the 
intensity: 

2 




X(r;m37r,t') = 


£=±1r=l 






(17) 


It is convenient to introduce the spin-density matrix 

= (18) 

r=l 


which represents the full information that is obtainable about X~. The diagonal elements of g 
are proportional to the partial-wave intensities and the off-diagonal entries to the interference 
terms. 

There are several effects that could lead to deviations from full coherence of the intermediate 
states, e.g. spin-flip and spin-non-flip processes or the excitation of baryon resonances. Also, 
performing the analysis over a large range of four-momentum transfer without taking into account 
the different t' dependences of the intermediate states may appear like incoherence (see Section . 
One way of including these incoherences is the introduction of the additional rank index r for the 
transition amplitudes, which is summed over incoherently [see Eq. The parameter Nr is 

called the rank of the spin-density matrix. 

The constraints due to parity conservation in the production process are directly taken into 
account by working in the so-called reflectivity basis, where positive and negative values for the 
spin projection M are combined to yield amplitudes characterized by M > 0 and an additional 
quantum number e = ±1, called reflectivity. This is achieved by replacing the D-function in the 
X~ —vr” two-body decay amplitude of Eq. 0 by 


'^-^ma(</'ty,'!?gj, 0) = c(M) 


-^(+m)a('^ty,^?gj,0) 


eP{-Y ^T)(^_^)^((/)ty,'!9gj,0) 


, (19) 


with e = ±1, M > 0, and the normalization factor 


c(M) = 


1/2 forM = 0, 
1 / \/2 otherwise. 


( 20 ) 
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The reflectivity is the eigenvalue of reflection through the X production plane. In the high-energy 
limit, £ corresponds to the naturality of the exchange in the scattering process, such that e = -\-l 
( —1) corresponds to natural spin-parity of the exchanged Reggeon, i.e. = 1~ or 2^ or 3“ 

... (unnatural spin-parity: = 0“ or 1"*“ or 2“ ...) transfer to the beam particle. Expressing 

the amplitudes in the reflectivity basis brings the spin-density matrix into a block-diagonal form 
with respect to e |29| . Hence states with different reflectivities, i.e. those produced by Regge- 
trajectories with different naturalities, do not interfere and are thus summed up incoherently 
[see Eq. (|17[)]. In general, the rank of the spin-density matrix may be different in the two 
reflectivity sectors, i.e. N^. 

Finally, we introduce the phase-space-normalized decay amplitudes as 








( 21 ) 


where dip^^r') is the differential three-body phase-space element. This normalizes the transition 
amplitudes via 

Ta%m3n, t') = Taimsn: t') \J ^ dip 3{T')\Wl{T'\m3^)\^ ( 22 ) 

with 


m 




(23) 


r=l 


such that the partial-wave intensities, which are the diagonal elements of the spin-density matrix 


in Eq. (22), are given in terms of number of events that would be observed in a perfect detector. 


The goal of the partial-wave analysis is to extract the unknown transition amplitudes ( 772 , 3 ^, t') 
from the data, because they contain information about the intermediate Svr resonances. Since the 
mass dependence of the transition amplitudes is unknown, the event sample is divided into 
bins much narrower than the width of typical hadronic resonances. Within each mass bin, the 
m 37 r dependence is assumed to be negligible, so that the 72^ only depend on t'. 

Also the t' dependence of the transition amplitudes is a priori unknown. In previous analyses it 
was often assumed that the m 3 ,r and t' dependences factorize and the t' dependence was modeled 
by real functions These functions were extracted from the analyzed data set by integrating 

over wide 7723 ^ ranges, often only for groups of waves. The COMPASS data, however, exhibit 
a complicated and significant correlation of the t' and 7723 ^ dependences (see Section]^, which 
renders this approach inapplicable. As it will be shown in Section |4.3| this is mainly due to 
different production processes (resonance production and nonresonant processes, like e.g. the 
Deck process |25]), which contribute with amplitudes that may have very different dependences 
on Therefore, the partial-wave decomposition is performed for each m^.,^ bin independently in 
different slices of t' (see Section 4.2 and Table®. Within a t' bin, the transition amplitude is 


assumed to be independent of t'. Taking out the explicit assumptions about the t' dependences 
by virtue of our large data set is an advantage compared to most previous analyses (e.g. 0). 

For a given bin in m^-,^ and t', the intensity has thus a simpler form as it depends only on the five 
phase-space variables r: 


m 


= E E 


£=±1r=l 




(24) 


with the transition amplitudes appearing as constants. Here, we introduced an additional 
incoherently added wave that is isotropic in r and from now on is referred to as flat wave. The 
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purpose of this wave is to absorb intensity of events with three uncorrelated pions in the hnal 
state, e.g. non-exclusive background. The flat wave is always part of the wave set, even if not 
mentioned explicitly. 


3.4 Maximum-Likelihood Method 

The transition amplitudes are determined for each bin in m^-,^ and t' by htting the model 
intensity T(r) of Eq. (24) to the measured r distribution. The ht is based on an extended 
likelihood function constructed from the probabilities to observe the N measured events with 
phase-space coordinates Ti'. 


C 


^ -N TT 

/d<^3('r)r/(r)X(r) ■ 

'-V-" '-V-' 

Poisson Probability for event i 

probability 


(25) 


Here, '^{t) is the detection efficiency and d(/93(r) the differential three-body phase-space element. 
The expected number of events N in the detector is given by the normalization integral 


N = 


dv93(r)r/(r)X(r). 


(26) 


By this integral, the detection efficiency is taken into account in the fit model, thereby avoiding 
the binning of the data, which would be impractical given the high dimensionality of the intensity 
distribution. 


Inserting Eq. (26) into Eq. (25) and dropping all constant terms as well as taking the logarithm, 
the expression reads 


N 


m 


2=1 *-£=±1 r=l 


+ 2flat 


tl r=l a 

EEEt;"’?’'' /<iV3(T)^(T)!P„'(T)4-r(T)-X«.. fd^ 3 (r)ri{T). (27) 

--LI ..._1 ^ J J 


£=±1 r=l a,b 


= 

— ^ab 




Here, the complex-valued integral matrix which is independent of the transition amplitudes, 
is calculated using Monte Carlo methods. The same is true for the real-valued integral /flat for 
the isotropic flat wave. 


In every individual bin, the transition amplitudes 72^ are determined by maximizing 


the likelihood function of Eq. (27), which allows the determination of the spin-density matrix 
elements 






(28) 


r=l 


Setting the detection efficiency ?/(r) 
number of events: 


1 in Eq. (26) gives the expected acceptance-corrected 
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Using the fact that the decay amplitudes normalized via Eq. (211 and that is 

hermitian, the expression can be rewritten as 


^ r.nrr — 


X] I I +2:flat f d^p3{r). 

e=±l as_. a<6s_1_ 


(30) 


CL a<h 

Intensities 


Overlaps 


From this equation, the interpretation of the spin-density matrix elements becomes obvious. The 
diagonal elements are the partial wave intensities, i.e. the expected acceptance-corrected 
number of events in wave a. The overlaps are the respective number of events that exhibit 


interference between waves a and b. Limiting the summation in Eq. (30) to a subset of partial 


waves yields the expected acceptance-corrected number of events in these waves including all 
interferences. Such sums will be denoted as coherent sums of partial waves in the following text. 

The procedure described in this section is referred to as mass-independent fit. It is worth stressing 


that fits in different kinematic bins are independent of each other. The fit model of Eq. (24) does 


not contain any assumptions about possible 37r resonances. They will be extracted in a second 
analysis step from the m 3 ,r dependence of the spin-density matrix. This so-called mass-dependent 
fit will be described in a forthcoming paper 


4 Partial-Wave Decomposition in Bins of and t' 


In principle, the partial-wave expansion in Eq. (24) includes an inhnite number of waves. In 


practice, the expansion series has to be truncated. This means that one has to define a wave set 
that describes the data sufficiently well, without introducing too many free parameters. 

Since the intermediate state X~ decays into a system of three charged pions, the G parity 
of X~ is —1 and the isospin 1=1, ignoring flavor-exotic states with / > 1. The number of 
possible partial waves is largely determined by the maximum allowed spin J of X~ , the maximum 
allowed orbital angular momentum L in the decay of the X~ to the isobar and the bachelor tt~, 
and the choice of the isobars. Since there are no known resonances in the flavor-exotic TT~Tr~ 
channel, we choose to include only 7 r~7r~^ isobars. We include [vrvrjs, p(770), /o(980), /2(1270), 
/o(1500), and /93(1690) as isobars into the ht model. This selection is based on the features 
observed in the Tr^vr"’' invariant mass spectrum in Figs. 5b and and on hndings of previous 
experiments uniiEiiiniEniES]. 

4.1 Isobar Parametrization 

In this section, we present the parametrizations of the mass-dependent amplitudes of the six 
isobars chosen above, which enter the analysis via Eqs. 0, (0, and and are summarized in 
Table |3l 

In most cases, the Tr~TT~^ isobar resonances are described using a relativistic Breit-Wigner 
amplitude m 

mo, Fq) = 

TTig — — I mo 1 (m) 

where mo and Tq are mass and width of the resonance. For a single two-body decay channel, the 
mass-dependent width r{m) is given by 

mo q F^{q) 


F{m) = Fo 


m qo Ff{qo) 


(32) 


By applying Eq. (32), we assume that the isobar decays predominantly into two pions and neglect 


other decay modes. Here, q{m) is the momentum of vr and vr'*' in the rest frame of the isobar 
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Table 3: Overview of the isobar parametrizations used in the partial-wave analysis. 


Isobar Formula 


Parameters 


[ttttJs 


p{77Q) 

/o(980) 


M solution from Ref. m 
(see Fig. 


101 


Eq. (31) with Eq. (40) 


Eq. (43) (see Ref. |l2]) 


/2(1270) 

/o(1500) 

P3(1690) 


Eq. 

Eq. 

Eq. 


(31) 


with Eq. 


(31) 


with Eq. 


(41) 


(32) 


(42) 


see text and Table 1 in Ref. m 

mo = 768.5 MeV/c^ 

To = 150.7 MeV/c2 

mo = 965 MeV / 
g-iT-K = 0.165 (GeV/c^)^ 

9kkI 9-k-k = 4.21 

mo = 1275.4 MeV/c^ 

To = 185.2 MeV/c2 

mo = 1507MeV/c^ 

To = 109MeV/c2 

mo = 1690MeV/c^ 

To = 190MeV/c2 


with mass m. At the nominal resonance mass, the breakup momentum is given by qo = q(mo). 
By Fi{q) we denote the Blatt-Weisskopf barrier factors |44| . which appear also in Eq. (|^ and 
take into account the centrifugal-barrier effect caused by the orbital angular momentum i in the 
isobar decayjj^ We use the parametrization of von Hippel and Quigg where 


FoiQ) 

Ffiq) 

Fiiq) 

Fiiq) 

Fiiq) 


1 , 

2z 

z + F 
13^2 

z^ + 3z + 9' 

277z^ 

z3 _|_ q^2 _|_ 45_2- _|_ 225 ’ 

12 746^"^ 

-p 10^3 -h 135^2 -h 1575z -h 11025’ 


(33) 

(34) 

(35) 

(36) 

(37) 


Fiiq) 

Fiiq) 


998 881^5 

-h 15^4 -h 315^3 -h 6300^2 -h 99 225z -h 893 025 ’ 

118 394 977^6 

76'+2i 75’ + 630? + 1^900^3 + 496 125^^ + g 823 27¥z + io¥o5l) 025 ’ 


(38) 

(39) 


Here, z = iq/qR)'^ with the range parameter qn = 202.4MeV/c that corresponds to an assumed 
strong interaction range of 1 fm 


For small breakup momenta ~ 0, the amplitude behaves like 


Fiiq) oc r 


The description of the p(770) isobar is slightly improved by modifying Eq. (32) as shown in 
Refs. |36l liT] : 

q Fjiq) 


rim) = To 


qo Ffiqo) 


(40) 


the decay of the isobar into two spinless particles, i is given by the spin of the isobar. 

Instead of the original normalization of the barrier factors such that Fi{q) —>■ 1 for q —>■ oo, von Hippel and 
Quigg modified the normalization in a way that Fi{q) = 1 for z = 1. 
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For the /?3(1690) isobar, a slightly modihed Breit-Wigner amplitude is used: 


^P3(1690)("^; "^0) -H)) — 


mo -To 


ttiq — 'm? — i 772-0 -Hd 


(41) 


Siuce the vr tt'*' decay mode of the /93(1690) is uot domiuaut, a coustaut total width is used. 

The most difficult sector is that of the scalar isobars with which cousists of several 

overlappiug /q resouauces. lu this aualysis, we cousider three iudepeudeut isobar amplitudes that 
have quite differeut properties. A broad compoueut with slow phase motiou, which we deuote by 
[vTTrj^, domiuates the mass spectrum from low to iutermediate two-piou masses. This compoueut 
iuterferes with the uarrow /o(980). lu elastic vrvr scatteriug, this iuterfereuce is destructive, so 
that the /o(980) appears as a prououuced dip. However, iu 7 r“ 7 r“ 7 r^ decays, the Trvr S-wave 
subsystem behaves differeutly. As will be showu iu Sectiou the relative phase betweeu the 
two compoueuts depends ou the quantum numbers of the Stt intermediate state and on its mass. 
In order to give the model the freedom to adjust the couplings of the various Svr states to the 
[ 7 r 7 r ]5 7 r and /o(980)7r decay modes separately, the broad tttt S'-wave component and the /o(980) 
are treated as independent isobars. Similarly, the /o(1500) is included using a Breit-Wigner 


amplitude [see Eq. (31)] with constant width 


r(m) = To. 


(42) 


The Breit-Wigner amplitude is not able to describe the /o(980) well as it peaks close to the KK 
threshold. Therefore, this isobar is described by a Flatte parametrization j48| that takes into 
account the coupling to the vrTT and KK decay channels: 




1 


777q — 772^ — i 




Qn-ir + 9kK 


(43) 


Here, is the two-body phase space for the two decay channels i = vrTT, KK with 

the respective breakup momenta qi(m), which become complex-valued below threshold. The 
values for the couplings and gj^j^ as well as that for the mass mo are given in Table as 
determined by the BES experiment from a partial-wave analysis of J/t/) decays into and 

(j)K-K+ @ 2 ]. 


The parametrization of the broad tttt S'-wave component is the most complicated one. It is based 
on the parametrization of the tttt S-wave from Ref. EH, which was extracted from vrTr elastic 
scattering data. We modify the so-called M solution (see Table 1 in Ref. [41j ) as suggested by the 
VES collaboration |49) . In order to remove the /o(980) from the amplitude, the parameters fl, 
/I) /i) C 33 , and C 22 as well as the diagonal elements of the M matrix in Eq. (3.20) of Ref. |41) are 

]). It has 


set to zero. Figure 10 shows the resulting [TTTrj^ amplitude (TIi of Eq. (3.15) in Ref. 
a broad intensity distribution that extends to two-pion masses of about 1.5 GeV/c^ accompanied 
by a slow phase motion. 


4.2 Fit Model 

When using the isobar model, we have in principle to take into account all observed 7 r“ 7 r'’“ 
correlations. In accordance with the vr^vr'’' invariant mass spectrum shown in Fig. |5b| and with 
analyses by previous experiments, we include [Trvrjs, p{770), /o(980), /2(1270), /o(1500), and 
P3(1690) as isobars into the ht model. Based on these six isobars, we have constructed a set of 
partial waves that consists of 88 waves in total, i.e. 80 waves with reflectivity e = -|-1, seven waves 
with e = — 1 and one noninterfering flat wave representing three uncorrelated pions. This wave set 
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Figure 10: Parametrization of the [ttttJs isobar amplitude based on the M solution described in Ref. im¬ 
panel (a) shows the intensity, (b) the phase, and (c) the corresponding Argand diagram. The open circles 
in the latter are evenly spaced in in 20MeV/c^ intervals. Arbitrary units are denoted by “a. u.”. 
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has been derived from a larger set of 128 waves, which includes mainly positive-reflectivity partial 
waves with spin J < 6, orbital angular momentum L < 6, and spin projection M = 0, 1, and 2. 
Omission of structureless waves with relative intensities below approximately 10“^ yields the 
88 partial waves that are used in this analysis and given in Table |A.l in Appendix 

The wave set includes waves with spin-exotic = 1 and 3 These waves have intensities 
significantly different from zero. They contribute 1.8% and 0.1 %, respectively, to the observed 
intensity. Removing the three 1 waves from the fit mode|i^ decreases the log-likelihood value, 
summed over the lit' bins, by more than 4000 units in the 37r mass range from 1.1 to 1.7 GeV/c^. 
If instead the two 3 waves are removed, the log-likelihood value, summed over the 11 t' bins, 
decreases by 200 units in the 37r mass range from 1.1 to 1.7GeV/c^. The spin-exotic waves will 
not be discussed any further in this paper. 


In the construction of the wave set, problems may arise when more than one isobar with the 
same quantum numbers and a broad overlap of their mass functions are used simultaneously, 
causing considerable overlap between the corresponding decay amplitudes. Such cases have to be 
treated with great care as the fit tends to become unstable. In our fit model, this applies to the 
O^'*' isobars. Here, the broad [vrvrjg, the narrow /o(980), and the /o(1500) do have considerable 
overlap. Because of the narrowness of the /o(980), the fit is able to separate it well from the 
broad [Trvrjs, as it is demonstrated in Section]^ In contrast, the inclusion of several waves with 
/o(1500) TT decay modes tends to destabilize the fit. Therefore, the 88- wave model includes only 
one /o(1500) vr wave. We decided to include the 0 O"'' /o(1500) vr S wave for 772 . 3 ^ > 1.7 GeV/c^ 

in order to study a potential signal for the decay 7r(1800) —)• /o(1500) -|- vr. The parametrizations 
used for the line shapes of the isobars are based on prior knowledge and were described in 
Section [ 4 .1[ The effect of isobars with uncertain line shapes may lead to spurious results and is 
addressed by systematic studies discussed in Section |4.6| and Appendix 
extended analysis method, which partly removes the model bias due to the isobar parametrizations. 
Results are presented in Section 


B.3 We also apply an 


The likelihood function to be maximized in the fit with the production amplitudes as free 
parameters is built according to Eq. (27|. Using such a large wave set to fit the three-pion 
system, we have to be concerned about stability of the results, which in turn may be influenced 
by correlations and cross talk of partial waves. In order to reduce such effects, different subsets of 
the 88 waves are used, which grow in size with increasing three-pion mass. High-spin waves and 
waves with heavy isobars are typically omitted from the wave set in the region of low msTr. This 
has two reasons: first, the intensity of such waves is expected to vanish at low msTr, and secondly, 
they would artificially contribute to ambiguities since the phase space at low masses appears 
to be too small to find a unique solution. A disadvantage of introducing the mass thresholds 
for particular waves are possible discontinuities induced in the mass dependence of other partial 
waves. Therefore, such thresholds have to be placed as low as possible. In our analysis, thresholds 
were applied to 27 of the 88 partial waves. The threshold values, which were carefully tuned in 

Appendix 


order to reduce artificial structures, are listed in Table A.l 


m 


The partial-wave analysis is performed independently in 100 equidistant msTr bins with a width of 
20MeV/c^, each of which is subdivided into eleven non-equidistant t' bins (see Table |^. The t' 
bins are chosen such that, except for the two highest t' bins, each bin contains approximately the 
same number of events. Within each of these 1100 bins, the transition amplitudes 72^{itistt, t' 


in 


Eq. (241 are assumed to be constant. Figure 11 illustrates the correlation of t' and rnsn, where 
the subdivision into bins of t' is indicated by horizontal lines. 


f^^This reduces the number of free parameters in the PWA fit by 6. 
I®'This reduces the number of free parameters in the PWA fit by 4. 
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OTj;! [GeV/c^] 


Figure 11: (color online) Correlation of the reduced four-momentum transfer squared t' and the invariant 
mass TO 3 ,r of the Stt system in the analyzed kinematic region. The partial-wave analysis is performed 
independently in 100 equidistant bins with a width of 20MeV/c^, each of which is subdivided into 
eleven non-equidistant t' bins. The latter are indicated by the dashed horizontal lines. The numerical 
values for the t' bins are listed in Table IH 


Table 4: Borders of the eleven non-equidistant t' bins, in which the partial-wave analysis 
(see Fig. 111. The intervals are chosen such that each bin contains approximately 4.6 x 10® 
the last range from 0.449 to 1.000 (GeV/c)^ is subdivided further into two bins. 


is performed 
events. Only 


Bin 


1 

2 3 

4 5 6 

t' [(GeV/c)2] 

0.100 

0.113 

0.127 

0.144 0.164 0.189 


Bin 


7 

8 9 

10 11 

t' [(GeV/c)2] 

0.220 

0.262 

0.326 

0.449 0.724 1.000 
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In the present analysis, we have limited ourselves to a rank-1 spin-density matrix for the waves 
in the posi tive- reflectivity sector, which is dominated by Pomeron exchange (see also discussion 
in Section 3.31. The smallest rank, = 1, is sufficient mainly because the analysis is 


performed in narrow bins of t', where the relative phases of the partial waves do not vary 
signihcantly. As part of the systematic studies, a ht with rank 2 was investigated, which shows 
enhanced artificial structures as well as increased instabilities (see Section 4.6 and Appendix |B.1[ ). 
For the negative-reflectivity waves, which can be produced by the exchange of various Reggeons 
[e.g. 5i(1235)], we use Nr^~ = 2. 


4.3 Selected Partial Waves with Spin Projections M = 0, 1, and 2 

In this section, we present the result of the fit in bins of m^-K and t' for 18 selected waves with 
positive reflectivity as listed in Table The waves are selected partly in view of the mass- 
dependent fit that will be described in a forthcoming paper, in which all resonance parameters 
determined by the fit will be presented |33]. This selection includes waves with spin projections 
M = 0, 1, and 2 that either have large intensities or exhibit clear signals of well-established 
resonances or even unexpected signals. In addition, we have selected large waves with [vTTrJg' and 
/o(980) isobars, which are related to the detailed study of the vr^vr'*' subsystem with = 0^"*“ 
presented in Section The amplitudes of 17 out of the 18 selected waves are found to be 
practically insensitive to systematic effects arising from the remaining waves, in particular to the 
truncation of the partial-wave expansion series in Eq. (24| (see Section 4.6 and Appendix [B|) . 


The intensity distributions of the remaining 69 waves are shown in the Supplemental Material in 
Appendix 

The total intensity of all partial waves is defined as the total number of acceptance-corrected 
events as given by Eq. (30). The relative intensity of a given partial wave, as listed in Table is 
defined as the ratio of its intensity integral over the analyzed mass range to the integral of the 
total intensity. This value is in general different from the contribution of a wave to the total 
intensity, owing to interference effects between the waves. Therefore, the relative intensities of 
all 88 partial waves add up to 105.3 % instead of 100 %. However, self-interference due to Bose 


symmetrization is included via Eq. (12). 


If not indicated otherwise, the wave intensities shown in the figures below are the sum of the 
intensities over the individual t' bins. They will be referred to as t'-summed intensities in the 
text that follows. The percent numbers given in the mass spectra are the relative intensities of 
the particular partial wave shown. In addition, we show for some waves the intensity distribution 
in individual t' bins. While mass and width of resonances do not depend on the production 
kinematics, coherent nonresonant contributions may vary in shape and phase with f. This may 
lead to significant t'-4ependent shifts of mass peaks. Examples of such effects are discussed below. 


As shown in Fig. |12[ waves with negative reflectivity, which correspond to unnatural-parity 
exchange processes, contribute only 2.2 % to the total intensity. The dominance of natural-parity 
exchange processes is expected at COMPASS energies because the Pomeron contribution is 
considered to be dominant. Therefore, we are only taking into account positive-reflectivity partial 
waves in the following. 

The incoherent isotropic flat wave turns out to contribute about 3.1% to the total observed 


intensity (see Fig. 13). This magnitude is roughly consistent with the background level that one 


expects from extrapolating the non-exclusive background component visible in Fig. 3b into the 
signal region. 


Figure 14 shows the t'-summed intensities of two major waves with spin projection M = 0, i.e. 
the l"*"^ O'*' p{770) TT S and 2 0^ /2(1270) tt S waves. Both exhibit clear peaks corresponding to 
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Table 5: Waves selected for presentation in this paper out of the much larger pool of 88 waves used in 
the mass-independent fit (see Table A.l| in Appendix [A|). The partial waves with [ttttJs and /o(980) isobars 
at the bottom of the table will be discussed in Section [6l The intensities are evaluated as a sum over the 
11 t' bins and are normalized to the total number of acceptance-corrected events. They do not include 
interference effects between the waves. 


tPC 


Isobaric decay Relative intensity [%] Shown in 


1++ 

0+ 

p{770) TT S 

32.7 

Figs. 

14a 

17a 

and 

17c 


1++ 

1+ 

p{770) TT S 

4.1 

Fig. 

.5a 




1++ 

0+ 

/2(1270) 7rP 

0.4 

Figs. 

16b 

20a 

and 

20c 


2++ 

1+ 

p{770)ttD 

7.7 

Figs. 

15b 

17b 

and 

17d 


2++ 

2+ 

p{770)ttD 

0.3 

Fig. 

.6a 




2++ 

1+ 

/2(1270) 7rP 

0.5 

Figs. 

15c 

20b 

and 

20d 


2“+ 

0+ 

p{770) TT F 

2.2 

Figs. 

21a 

21c 

and 



2“+ 

0+ 

/2(1270) 7r5 

6.7 

Figs. 

14b 

18a 

18c 

anc 

22a 

2-+ 

1+ 

/2(1270) 7r5 

0.9 

Figs. 

15d 

and 

22c 


2-+ 

0+ 

/2(1270) 7rL» 

0.9 

Fig. 

22b 




4++ 

1+ 

p{770)ttG 

0.8 

Figs. 

15e 

18b 

and 

18d 


4++ 

1+ 

/2(1270) 7rF 

0.2 

Figs. 

15f 

21b 

and 

21d 


0-+ 

0+ 

[vrvrjs TT S 

8.0 

Figs. 

24a 

24b 

and 

25b 


0-+ 

0+ 

/o(980) 7r5 

2.4 

Fig. 

25a 




1++ 

0+ 

[TTVrj^ TT P 

4.1 

Fig. 

25f 




1++ 

0+ 

/o(980)7rP 

0.3 

Fig. 

25e 




2-+ 

0+ 

[7r7r]5 TT D 

3.0 

Fig. 

25d 




2“+ 

0+ 

/o(980)7rD 

0.6 

Fig. 






Intensity Sum 


75.8 
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Figure 12: The i'-summed intensity of the coherent sum of all negative-reflectivity waves (a) and, for 
comparison, together with the total intensity of all partial waves (b). 
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Figure 13: The t'-summed intensity of the flat wave (a) and, for comparison, together with the total 
intensity of all partial waves (b). 
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xio^ 1+Wp(770);r5 




Figure 14: The i'-summed intensities of major waves with spin projection M = 0 showing in (a) the 
ai(1260) and in |(b)| the 7r2(1670). The shaded regions indicate the mass intervals that are integrated over 
to generate the t' spectra (see Fig. 32a I. 


the ai(1260) and the 7r2(1670) resonances. 

Selecting spin projection M = 1, we have access to the I"*" p{770) vr S, 2++ 1+ p{770) tt D, 2++ 
/2(1270) vr P, 2 ^ f2{1270) tt S, as well as to the 4'^^ p{770) tt G and I"''/2(1270) 



agreement with our result for the 2'^~^ wave in the r/vr final state, which is dominated by the 
02 ( 1320 ) [50]. In the analyzed range of 0.1 <t'< 1.0 (GeV/c)^, the observed suppression is twice 
as large as the suppression of M = 1 versus M = 0 waves. 


Nonresonant and resonant contributions are expected to follow different production paths with 
possibly different dependences on t'. In order to investigate possible nonresonant contributions, 
we show in Figs. 17 and 18 the intensities of four selected waves for two intervals of t', i.e. 
0.100 < t' < 0.113 (GeV/c)^ and 0.449 < t' < 0.724 (GeV/c)^, which represent regions of low 
and high f in this analysis. When comparing these two regions, the shapes of the 02(1320) and 
04 ( 2040 ) resonances in the 2'*“'’'1^ p(770) vr P and d"*"*'I"*"/?(770) tt G waves, respectively, stay 
largely unaltered. In contrast, we observe that the peak in the 1+'’' O"*" /?(770) tt S wave, which 
presumably contains the oi(1260), significantly shifts towards higher masses with increasing 
t'. A similar but less strong effect is observed for the 7r2(1670) peak in the 2 ^ 0"*“ /2(1270) tt S 
wave. This shows that the peak structures in the latter two partial waves are not only due to 
ordinary resonances but are distorted by nonresonant contributions. The Deck process proposed 
in Ref. |25| and illustrated in Fig. 19 may provide an explanation for the t'-dependent nonresonant 
contributions observed in the l"*"^ and 2 ^ waves. The t' dependence of the shape of the I"’""*' O'*’ 
p(770) TT S mass spectrum was already observed by the ACCMOR collaboration |141 ITH] and our 
results confirm their findings. 


We show in Figs. 20 and 21 the same t' regions for the small-intensity waves I’*"*' 0’*’ /2(1270) ttP, 
2’*'+ I"*" 72 ( 1270 ) TT P, 2 *■ O’*" p{770) vr F, and d'*"*' I"*" /2(1270) tt F. All waves show a pronounced 
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Figure 15: The ^'-summed intensity of waves with spin projection M = 1 showing in (a) the ai(1260), in 


1(b) I and [(c) I the 02(1320), in |(d)| the 7r2(1670), and in (e) and |(f)| the 04(2040). The shaded regions indicate 
the mass intervals that are integrated over to generate the t' spectra (see Figs. 32b 33a 33b and 341. 
































30 


> 

u 


o 

(N 


40 


30 


20 


10 


xlO' 

_0.3% 


2++2"^p(770) TtD 
0.100 < f'< 1.000 (GeV/c)2| 


(a) 


3.5 


1.5 


2 2.5 

[GeV/c^] 


o 

> 

u 


o 

<N 


C 

(U 


xio" 


l++0+/^(1270) TIP 


10 


0.4% 0.100 < f' < 1.000 (GeV/cr 




.V 




(b) 


3.5 


1.5 


2 2.5 

mjjj. [GeV/c^] 


Figure 16: Panel (a): t'-summed intensity of 2++ 2+ p(770) tv D wave with spin projection M = 2 and 
the 02 ( 1320 ) peak. Panel (b): t'-summed intensity of 1++O’*"/2(1270) tt P wave with spin projection 
M = 0. In this wave, the mass regions below 1.5GeV/c^ and above 2.0GeV/c^ (shown by gray points) 


are sensitive to the truncation of the partial-wave expansion series (see Section 4.6 1 . In both waves, the 


shaded regions indicate the mass intervals that are integrated over to generate the f spectra (see Figs. 33c 


and 33c I. 


dependence of the mass spectrum on t'. In contrast to the p{770) tt S wave, the Oi(1260) 
cannot be clearly identihed in the /2(1270) tt P wave. Instead, the latter wave shows a broad 
enhancement around 1.8GeV/c^ (see also Fig. 16b). In the /2(1270) vr P wave, the 02(1320) 


exhibits a high-mass shoulder, which is particularly pronounced at large values of t', although it 
is clearly identihable also at low Such a high-mass shoulder also becomes prominent for the 
7r2(1670) in the p(770)7rP wave, for which the spectrum exhibits a richer structure than for the 
2 ^ O'*" 72 ( 1270 ) vr S wave. 

The selective effect of the orbital angular momentum L in the decay is clearly demonstrated in 
Figs. 22a and 22b which show the 2 ^ O"*" /2(1270) vr waves with L = 0 and L = 2. The 7r2(1670) 


dominates the S'-wave, while the 7r2(1880) favors the D-wave. The 7r2(1880) is considerably lighter 
than the expected radial excitation of the 7r2(1670) ground state and has been rated as a viable 
hybrid-meson candidate PEI]. However, a dominance of S over D-wave / 2 ( 1270 ) vr decay modes 
was predicted for hybrid mesons by model calculations |52[ I53| . This is at variance with the 
present observation of a prevailing D-wave decay of the 7r2(1880). The existence of the 7r2(1880) 
was questioned by Ref. |54| . which explains it as an interference of the 7r2(1670) ground state 
with the nonresonant Deck process causing an apparent shift of the peak position. This might be 
counterargued by the observation of two peaks in the 2 O"*" p{770) tt F wave (see Figs 


21a 


21c 


and 22d). Studying additional decay modes and the f dependence of the mass spectra should 


resolve this issue. 


4.4 Partial Waves with /o(980) and broad tttt S^-wave Isobars 

As discussed in the previous section, the shape of the peak in the ai(1260) region in the l"*"^ O'*" 
p(770) TV S wave and that in the 7r2(1670) region in the 2 ^ waves change as a function of t' (see 


Figs. 17a, 17c 18a and 18c|. A possible explanation for this behavior is the Deck process. We 


have therefore investigated partial waves that are expected to have small contributions from the 
Deck process. Owing to the nature and small width of the /o(980), this is in particular true for 
/o(980) TT partial waves. Only a few meson resonances have been observed to decay via /o(980), 
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Figure 17: The intensities of the I++O'*'p(770) tt S' and 2++1+p(770) tt Z3 waves in two different t' 
regions. Upper row: low t'; Lower row: high t'. The shaded regions indicate the mass intervals that are 
integrated over to generate the t' spectra. 
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Figure 18: Same as Fig. but for the 2 0+ /2(1270) tt S and 4++ 1+ p(770) tt G waves. 
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Figure 19: Example for a nonresonant production process for the Stt final state as proposed by Deck |25| . 
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Figure 21: Same as Fig. but for the 2+0+ p{770) tt F and 4++ 1+ /2(1270) tt F waves. 






















35 


xlO® 


2 0+/^(1270) TtS 


xlO" 


2~W/^(1270) tzD 


> 


o 

(N 


C 

u 


0.3 


0.2 


0.1 


6.7% 0.100 <f' 

< 1.000 (GeV/c)^ 


0.9% 0.100 <t’< 1.000 (GeV/c)^ 



> 

- 


- 

. 

§ 40 

- 


- 


w 



- 




* 

- 

- 








. 


* 

C 

” 20 

* 







1 i 1 



(a) 


}.5 


xlO" 


1.5 2 2.5 

[GeV/c^] 


2-nV2(1270) ns 


(b) 


3.5 


xlO" 


1.5 


2.5 


[GsV/c^] 
2" W p(770) n F 


o 

> 


o 

(N 


C 

u 


40 


20 


.0.9% 0.100 <f' 

< 1.000 (GeV/c)^ 

^ 80 

2.2% 0.100 <f' 

< 1.000 (GeV/c)^ 




> 








1 

- 







o 

07) 

- 

* 








♦ 



+ 






+*♦ ♦ 


• 


• 

(/3 

§ 40 

- 




,* 


* 

c 





•* 


• 

20 

- «♦ 



**.* 








****•. 

*“*’**'** 




. 




5 1 1.5 

2 2.5 5 

5 1 1.5 



1 2 


[GeV/c^ 

(d) 



'”3;r 

[GeV/c2] 


(c) 

Figure 22: The i'-summed intensities of 2 ^ M+ /2(1270) n L waves in panels (a), (b), and (c) and of 
2 O’*" p(770) TT F wave in panel (d). The 7r2(1670) dominates the waves shown in (a) and (c) the 7r2(1880) 


the one shown in |(b)[ Both resonances appear in the wave shown in |(d)[ Upper row: comparison of same 
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Figure 24: Intensity of the 0 “*■ 0+ [ttttJs tt S wave in two different t' regions, (a) low t'; (b) high f'. The 
shaded regions indicate the mass intervals that are integrated over to generate the t' spectra. 


such as 7r(1800) ^ /o(980)7r, (/.(1020) ^ /o(980)7, (/>(2170) ^ /o(980)(/>, and 77(1405) ^ fo{980)r], 
where the latter is a sub-threshold decay. Among these, the 7r(1800) is the only isovector state 
and thus accessible in the Stt final state. Because of its small width, the /o(980) accounts for 
only a small fraction of the full vrvr S'-wave. It is easily separated from the broad vrvr S'-wave 
structure, which is shown in Fig. Compared to the positive-reflectivity waves containing the 
p{770) isobar, those with the /o(980) are suppressed by a factor of approximately 20. Figure 23 
shows the t'-summed intensity of the coherent sum of all partial waves with an /o(980) isobar 
and positive reflectivity, which amounts to a relative intensity of 3.3 %. 

The intensity distribution of the 0 O'*" /o(980) vr S wave is dominated by the 7r(1800) peak (see 

Fig. 25a|. The more complicated mass spectrum of the 2 O"*" /o(980) vr Z3 wave is shown in 


Fig. 25c This wave, which should contain signals of 7r2(1670) and 7r2(1880), is characterized by 
pronounced destructive interference around m^.,^ = 1.8GeV/c^. 


The I"*"*" O"*" /o(980) vr P wave is shown in Fig. 25e It exhibits the new axial-vector meson ai(1420) 
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Figure 25: The t'-summed intensity for selected waves with tttt S'-wave isobars. Left column: waves with 
the narrow /o(980) isobar; Right column: waves with the broad isobar. The = 0 waves in 
the top row show the 7r(1800). The structure at 1.2 GeV/c^ is probably mainly of nonresonant origin. The 
2 intensities in the center row exhibit a complicated destructive interference pattern around 1.8 GeV j<?. 
The bottom row shows an enhancement in the region of the ai(1260) in the 1“'""'" 0"*" [ttttJs ttP wave and a 
new state, the ai(1420), in 1++ O'*' /o(980) ttP. The shaded regions indicate the mass intervals that are 


integrated over to generate the t' spectra (see Figs. 35 and 361. 
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The resonance features of this signal were presented in Ref. j55) . It should be noted that the 
intensity of this wave corresponds to only about 0.3 % of the total number of events. Since the 


signal is very small, we conducted several systematic studies that will be discussed in Section 4.6 


The same partial waves discussed above are shown in the right column of Fig. 25 for the broad 


component of the vrTr S'-wave as isobar, which is parametrized as described in Section 4.1 and 
denoted by [vrvrjs'. The intensity spectrum of the 0 O"*" [7r7r]5 7r5 wave (see Fig. 25b) exhibits 

two pronounced maxima and differs considerably from that of the corresponding /o(980) vr wave 
in Fig. 25a, The maximum at 1.8GeV/c^ corresponds to the 7r(1800), but we also observe a 
broad structure around 1.2GeV/c^, which could contain the 7r(1300). As it will be discussed in 


Section |5.2[ the latter structure has a very distinct dependence on t' with a minimum around 
0.35(GeV/c)2. 


t' 


Figure 


24 


shows as an example the intensity spectrum in two t' bins. At high 
t', the 7r(1800) peak nearly vanishes and the structure around 1.2GeV/c^ is shifted towards lower 
TTT-sjr. A more detailed analysis discussed in Section 6.3 indicates that in addition to interference 
effects with the 7r(1800) also nonresonant processes seem to contribute to the 1.2GeV/c^ mass 
region. 

Another resonance that was observed to couple to the [vrvrjs' isobar is the 7r2(1670) [9], whose 


main decay mode into /2(1270) tt is discussed above in Section 4.3 The mass spectrum of the 
2 ^ 0^ [7r7r]5 7rZl wave is shown in Fig. 25d and exhibits marked destructive interference effects 


at masses around 1.8GeV/c^, similar to the ones observed in the corresponding wave with the 


/o(980) decay mode in Fig. 25c 


The 1 O"*" [vrvrjs ttP wave is even more difficult to interpret (see Fig. 25f). A significant signal 

is observed in the region of oi(1260). However, as it will be shown in Section 5.2 and Section]^ 
this structure exhibits a strong t' dependence, which is a signature for signihcant nonresonant 
contributions. 

Comparing the /o(980) 7r and [Trvrj^Tr decay modes, the latter are obviously more affected by 
nonresonant contributions. We will discuss the sector of partial waves with TTTr S-wave isobars 
again in Section in the context of an extended analysis. 


4.5 Comparison of Fit Result and Real Data 

In order to estimate the goodness of the mass-independent ht, three-pion phase-space Monte 
Carlo events, which were processed through the detector simulation and reconstruction chain 
and satisfied the selection criteria, were weighted with the intensity distribution of the fit model 


[see Eq. (24)]. For a good fit, distributions obtained from these weighted Monte Carlo events are 
expected to approximate the real data. 

For fixed values of msTr and t\ the phase space of the three final-state particles is five-dimensional. 
Therefore, we can show only projections in certain kinematic regions. For the comparison we use 
the same hve kinematic variables that also enter in the decay amplitudes (see Section 3.2), i.e. 


cos??GJ cind (/)ty of the isobar in the Gottfried-Jackson frame, the isobar mass mg = and 

cos dfiF and ij^hf of the tt~ in the helicity frame. 


Figures 26 to 29 show as examples the distributions of the kinematic variables in various regions 


of msjT and t'. These kinematic regions contain different resonant and nonresonant contributions 
leading to different shapes of the angular distributions and the isobar mass spectrum. 

In general, the agreement between the weighted Monte Carlo and the real-data events is very good, 
in particular at intermediate 37r masses. At larger m^.^, we observe small deviations concerning 
the description of the /o(980) and /2(1270) isobars (see Fig. 27b) as well as small localized 


differences in the angular distribution in the Gottfried-Jackson frame (see Figs. 27e and 29e) 




























39 


0.127 < f'< 0.144 (GeV/c)2 
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Figure 26: (color online) Distributions of the five phase-space variables used to calculate the decay 
amplitudes shown for different Stt mass bins in the region 0.127 <t' < 0.144 (GeV/c)^. Each distribution 
is shown for real data (blue points) and for weighted Monte Carlo events (red bands), which are generated 
according to the fit result. Each distribution is normalized to its maximum deviation from its average 
y value. Along the ordinate, the average y values for the distributions (indicated by gray lines) are shifted 
equidistantly with respect to one another. The Stt mass ranges from 1.0 to 2.4GeV/c^ and is subdivided 
into 40MeV/c^ wide bins. The central values of selected Stt mass bins are given as labels in the cosdcj 
distribution. 
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0.127 < f'< 0.144 (GeV/c)^ 1.6 < < 1.8 GeV/c^ 




Figure 27: (color online) Comparison of kinematic distributions of weighted Monte Carlo events, generated 
according to the fit result, with the corresponding real-data distributions in the low-t' region 0.127 < t' 
< 0.144 (GeV/c)^. Panel (a) shows the acceptance-corrected Stt invariant mass distribution. The other 
panels show kinematic distributions in the mass interval 1.6 < < 1.8GeV/c^ around the 7r2(1670), 

which is indicated by vertical red lines in (a), (b) invariant mass spectrum of the 7r“7r+ subsystem; 
(c) distribution of the Gottfried-Jackson angles for real data; (e) ratio of the real-data distribution in (c) 
and that of the weighted Monte Carlo. Panels (d) and (f) show the respective distributions for the helicity 
angles. Note that (b), (c), and (d) have two entries per event. 
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Figure 29: (color online) Same as Fig. 27 but for the high-t' region 0.449 < t' < 0.724 (GeV/c)^ and 
the mass slice 1.2 < < 1.4GeV/c^ around the 02(1320). 
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4.6 Systematic Studies 


Given the high precision of the data, statistical uncertainties are negligibly small in most cases, 
i.e. systematic uncertainties are dominant. We have performed several tests to assess the stability 
of the result of the mass-independent fit. Here, we summarize the findings of these studies; more 
details can be found in Appendix [B| 

Possible effects from processes, in which the target proton is excited, are expected to be negligible. 
Due to Pomeron dominance, target excitations will be mostly N*. The recoil-proton trigger and 


the momentum-conservation criterion applied in the event selection (see Sections 2.2 and 2.3) 
suppress such events on average by about an order of magnitude. The remaining contributions 
consist predominantly of low-mass N* produced at large t'. In diffractive reactions, target and 
beam vertex factorize, so that these events are expected to have only little effect on the production 
of the three-pion final state. As we assume for these events the proton mass for the mass of the 
recoiling particle, the calculated values of t' would be slightly shifted by values comparable to the 
t' resolution. The same is true for the reconstructed beam energy illbeam- 

In order to reduce the probability for the fit to converge to a local maximum, the likelihood fit 
is repeated in each kinematic bin in msjr and t' 30 times with random starting values for the 
transition amplitudes in Eq. (27). From these 30 fits, the one with the highest likelihood is 


selected in order to generate the results presented in Sections 4.3 and |4.4| In the 37r mass range 
above about 1 GeV/(?, the fits reliably yield a single solution. Only a few fits are trapped in local 
maxima with significantly lower likelihood. In contrast, for mass bins below about IGeV/c^, the 
fits find multiple local maxima that deviate from each other only by a few units of log-likelihood. 
We attribute this behavior to the fact that, due to the smaller phase-space volume at low msTr, 
mainly the low-mass tails of the isobars contribute to the decay amplitudes. Therefore, it is harder 
to distinguish partial waves with different isobars. Since we do not expect any 37r resonances 
below 1 GeV/c^, no efforts were made to resolve these ambiguous solutions. 


We have studied how the truncation of the partial-wave expansion series in Eq. (24) influences 
the intensities of the 18 partial waves discussed in this paper (see Table by comparing to a 
fit with a reduced set of only 53 partial waves |56) . Except for one wave, the intensities exhibit 


only relatively small changes, which typically affect the high-mass regions. The intensity of the 
0+ /2(1270) vrP wave changes significantly in the mass regions above 2.0GeV/c^ and below 
1.5 GeV/c^, the latter of which is attributed to model leakage from the /?(770) vr S-wave decay 
of the ai(1260). However, the region around the enhancement at 1.8GeV/c^ (marked by the 


shaded region in Fig. 16b) is only slightly affected 


As mentioned in Section |3.2[ we do not apply relativistic corrections to the decay amplitudes in 
the partial-wave analysis. First studies show that the effect on the shapes of the selected 18 waves 
is small. 


In order to study the effect of the rank of the spin-density matrix, the data are fit with a rank-2 
spin-density matrix instead of the rank-1 for the positive-reflectivity sector. The most striking 
feature of the rank-2 fit, which has nearly twice the number of free parameters, is that the 
flat wave practically disappears. In addition, intensity is shifted from the negative into the 
positive-reflectivity sector. However, the fit shows substantial instabilities and artificial structures 
in the msTr region between 1.0 and 1.3GeV/c^ in some partial waves. We therefore conclude that 
the rank-1 fit offers a better description of the data using significantly less parameters. For more 
details see Appendix |B.l 

Omitting the waves with negative reflectivity from the wave set feeds intensity mostly into the 
flat wave, but causes little change of positive-reflectivity waves (see Appendix |B.2[). This means 
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that in the given range of four-momentum transfer, the positive and negative-reflectivity sectors 
do not mutually influence each other and are well-separated by the fit as opposed to the case of 
very small t' < 10“^ (GeV/c)^ |7]. 


The isobar parametrizations (see Section 4.1) are an important input for the PWA model. The 
p(770) is the dominant isobar. For most of the waves in Table the intensity distribution is not 
sensitive to the details of the p(770) parametrization or to small changes of the used parameter 
values. In contrast, the region around the Oi(1260) mass in the O'*' [7r7r]5 7rP wave and the 
02 ( 1320 ) signal in the 2'^^ 1+ /2(1270) vrP wave change significantly (see Appendix B.3|. Both 
seem to be contaminated by model leakage due to the imperfect description of the /9(770). The 
dependence of the PWA result on the parametrization of the /o(980) and that of the broad vrTr 
S'-wave component is studied as well. Using a simple S-wave Breit-Wigner amplitude [Eq. ( |31[ ) 
with Eq. ( [ 42 ) )] instead of the Flatte parametrization [Eq. ( |43| )] for the /o(980) reduces the height 
of the intensity peaks of the resonances decaying into /o(980)7r by about a factor of two (see 
Appendix |B. 3 j ). However, the shapes of the peaks in these partial waves remain unchanged. In the 
mass region above 1.3GeV/c^, the fit with the /o(980) Flatte parametrization has a significantly 
higher likelihood than the one with the S'-wave Breit-Wigner amplitude. This indicates that the 
data are better described by the Flatte parametrization. 


The influence of the parametrization used for the broad component of the tttt S-wave on the 
PWA result is studied by performing a fit with an alternative description of the mass-dependence 
of the isobar amplitude. 


Instead of the modified M solution from Ref. |41] (see Section 4.1), the 
Ki solution from Ref. m with the /o(980) pole subtracted using a simple S-wave Breit-Wigner 
amplitude is used. This parametrization was originally used by the VES experiment j20) . The 
Breit-Wigner for the /o(980) is similar to the one employed in the /o(980) study described above. 
In order to be consistent, the Breit-Wigner amplitude is also used for all waves with the /o(980) 
isobar. The observed variations in the fit result are small. 


Performing the PWA on a data sample without the GEDAR, RIGH, and central-production vetos 
described in Section 2.3 shows that the result is not very sensitive to backgrounds from kaon 
diffraction, kaon pairs in the final state, and central-production reactions (see Appendix |B.4| ). 
The partial-wave intensities scale approximately with the number of events, only the relative 
intensity of the flat wave increases. It is very unlikely that the peak around 77137 ^ = I.4GeV/c^ in 
the I"*"*' 0^ /o(980) tt P wave is caused by kaon-induced reactions or that it stems from kaonic 
final states misinterpreted as pionic ones. The GEDAR and RIGH vetos applied in the event 
selection reduce such contaminations considerably. Eurther studies show that the signal is not 
correlated with these cuts. 


In summary, the PWA fits converge reliably for flvr masses above about I GeV jc?. The shapes of 
clear resonance peaks are stable with respect to changes of the PWA model. However, in some 
cases the height of the intensity peaks is sensitive to the isobar parametrization. This issue will 
mostly be resolved for the Trvr S-wave isobars by applying a method introduced in Section § by 
which the isobar amplitudes are extracted from the data. 


5 t' Dependences 


Figures 30a and 30b illustrate how the shape of the measured three-pion invariant mass spectrum 
changes with t', while Figs. 30c and 30d show the change of the measured t' distribution with 
37r mass. It is apparent that the t' spectrum strongly depends on This observation has 

motivated us to perform the PWA in bins of t' instead of weighting the partial waves with 
t'-dependent model functions, which has been the conventional approach. 
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Figure 30: The t' dependence of the measured Stt invariant mass spectrum and vice versa. 
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Figure 31: (color online) Result of a fit to the t' dependence of events from the diffractive-dissociation 
reaction 7r“ + p —>■ 7r“7r“7r+ + p, as measured by COMPASS (filled circles) and by the ACCMOR 
collaboration [16] (open circles). For each 7713 ^ bin, the t' spectrum was fit using a double-exponential 
model [see Eq. (441]. Panel (a) shows the mass dependence of the two slope parameters bi (upper points) 
and &2 (lower points). Note the extended mass range of the present measurement as compared to the 
ACCMOR data. Panel (b) shows the ratio I 1 /I 2 of the integrated exponential contributions, see Eq. (451. 


In this chapter, we will elaborate on the details of the observed t' distributions. We will start with 
global spectra, from which we determine the t' dependence as a function of Such an analysis 
comes closest to the traditional description of high-energy reactions in terms of Regge exchange. 
Using the results from the mass-independent fit outlined in Section]^ we can in addition separate 
the contributions of various partial waves to the t' spectrum. Comparing different Svr mass regions 
that are either dominated by well-established resonances or by nonresonant contributions, various 
patterns become apparent. 


5.1 Overall t' Dependence 

The first extensive study of the Stt mass dependence of the t' spectrum was performed by the 
ACCMOR collaboration |13l IT 6 ] . They investigated the reaction 7 r“ -|-p —)• 7 r“ 7 r“ 7 r'*' -|-p at 
63 and 94 GeV/c incoming pion momentum and determined the t' dependence as a function of 
msjr in the range t' < 1.0(GeV/c)^. The t' dependence was parametrized for each 50MeV/c^ 
wide 37 r mass bin by two exponentials: 


diV 

W 




(44) 


with real-valued parameters A,. The ACCMOR collaboration observed that the two slope 
parameters 61^2 are different at small values of 7713 ^ and that they vary significantly up to values 
of m 37 r of about 1.2GeV/c^. This marks the onset of resonance production. Beyond this mass 
value, the slope values of 61 ~ 12 (GeV/c)“^ and 62 ~ 5 (GeV/c)“^ stay almost constant (see 


open circles in Fig. 31a). 


We perform the same study on the present data in a wider Sir mass range and using finer 772 . 3 ^ bins 
of 20 MeV/width. In each mass bin, the acceptance-corrected f spectrum, which is obtained 


from the mass-independent fit in 11 t' bins, is fit by Eq. (44). The result is shown as filled circles 
in Fig. |31a 


The general pattern and also the absolute values for the slope parameters agree 
nicely with the ACCMOR results. We observe a strong dependence of both slope parameters on 
















47 


msTT- In the region m3,r < l-OGeV/c^ below the resonances, only a few partial waves contribute 
significantly to the spectrum, i.e. the 1++ O’*" p(770) vr S (57.1 %), 0 ^ O’*" [vrvrjs vr S (12.8 %), l^^ 
p{770)ttS (8.3%), 0+[vrTr ]5 TT P (4.2%), and 0 ^ 0'^ p{770) tt P (4.1%) waves. In this 

region, the parameter bi representing the steeper component shows a rapid drop with increasing 
three-pion mass. The parameter 62 representing the shallower component exhibits less variation. 
Its mass dependence shows a dip by 2 (GeV/c)“^ at around = l.OGeV/c^. Approaching 
the mass region of ai(1260) and 02(1320), above approximately 1.3GeV/c^, the msTr dependence 
of bi and 62 changes abruptly: bi drops much slower, decreasing from about 12(GeV/c)“^ at 
1.3 GeV/c^ to 8 (GeV/c)“^ at 2.5 GeV/c^, whereas 62 stays nearly constant at about 4 (GeV/c)“^ 
over the same mass range. 


Figure [Mb shows the ratio of the contributions 


h = A, 


dt' 


— 61 t' 


and I 2 = A 2 


dt' 


=-b2 t' 


(45) 


of the two exponentials from Eq. (44), integrated from = 0.1 (GeV/c)^ to = 1.0(GeV/c)^ 


As observed for the slope parameters, the regions below and above the resonances show very 
different behavior. Below the resonance region, the component with the steep slope bi dominates 
and its contribution reaches a maximum at 37r masses of approximately 1.0 GeV/ (P‘. From there, it 
drops quickly with a shallow minimum around m3,r = 1.3GeV/c^. This dip is presumably caused 
by the onset of the 2~^~^ waves with M = 1. Above about 1.3GeV/c^, the relative contributions 
of the two exponentials only depend weakly on 772,3^ with almost equal relative weights for the 
two terms. The AGGMOR data show a qualitatively similar behavior. The agreement with the 
present data is, however, not as good as that observed for the slopes. 

To our knowledge, the complicated mass dependence of the t' spectra described above is not 
well understood. In the region around 1.3GeV/c^, nonresonant processes are known to play 
an important role. Most available calculations describe these processes as the dissociation 
of the beam pion into the isobar and the bachelor 7r“, followed by diffractive scattering 
of one of the beam fragments (typically the vr”) off the target proton (see Fig. 19). These 


calculations focus mainly on the 37r mass dependence and are based on vrvr and vrp elastic¬ 
scattering data I25i Ea EH ca Eg. The more elaborate three-component Deck model |59[ [6OI IGT] 


describes the reaction vr” -|-p —>• vr” -|-Precoil by including ^ as well as vr exchanges in addition to 
direct production of vr” via Pomeron exchange. In general, such nonresonant processes exhibit 
a dependence on t! that is different from that of resonant production. Interferences between 
resonant and nonresonant processes may in addition modify the t! spectra. The three-component 
Deck model describes the correlations between the ^vr invariant mass, the slope of the t' spectrum, 
and the cos v^gj distribution in detail and predicts the existence of interference minima in the t' 
spectra. As it will be shown in the following section, the t' spectra of some partial waves exhibit 
such kind of minima in certain 3vr mass regions. 

5.2 t' Dependences of Individual Partial Waves 

Using the partial-wave decomposition of the mass spectrum from the mass-independent fit as 


presented in Section 4.3 we can now study the t' dependence of the intensity of individual partial 


waves in different mass regions. The selected mass regions are indicated by shaded bands in the 
spectral distributions shown in Section The corresponding t' spectra are obtained by integrating 
the partial wave intensities over those mass regions. The integrated intensities are presented using 
horizontal bars, the lengths of which represent the widths of the given t' bins. Blue horizontal 
lines represent the central values. The height of the gray horizontal bars corresponds to the 
statistical uncertainty of the intensity. 
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Figure 32: (color online) The t' dependence of the intensity of the M+ p(770) tt S waves with spin 
projections M = 0 (a) and M = 1 (b), integrated over the mass region around the ai(1260) as indicated 


by the shaded regions in Figs. 14a and 15a The solid red curve in (a) represents a double-exponential fit 


using Eq. (461, the one in (b) a single-exponential fit with parameter A 2 = 0. In (a), the two exponential 


components are shown by the dotted lines. In (b), the fitted t' range is indicated by the solid curve; the 
extrapolation is shown as dashed curve. See text for details. 
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Figure 33: (color online) The t' dependence of the intensity of three 2++ waves with different isobars 
and different spin-projections M, integrated over the mass region around the 02(1320) as indicated by the 
shaded regions in Figs. |15bl |15c[ and |16a| The solid red curves represent double-exponential fits using 


Eq. (46); the fit components are shown as dotted curves. See text for details. 


We compare waves with the same isobars and angular momentum in the decay but with different 


spin projections M. Figure 32 gives an example for the M"*" /3(770) tt S intensities with 
M = 0 and 1 integrated over the range 1.1 < msTr < 1.3GeV/c^, which covers part of the 
ai(1260). Figures 33a and 33c show the analog comparison for 2'’”*“ M"*" /9(770) tt D waves with 
M = 1 and 2. 


Alternatively, waves with same quantum numbers but different isobars can be compared, again 


keeping the mass interval fixed. Figure 33 shows this for the 2'^~^ I"*" waves with the p(770) 
isobar and L = 2 along with waves with the /2(1270) isobar and L = 1, both in the mass region 
1.2 < m 3 ,r < 1.4GeV/c^ around the 02(1320). A comparison of the t' spectra of the 4+'’' I"*" 
p(770) ttG and /2(1270) tt F waves for 1.86 < < 2.06GeV/c^ is given by Fig. 


34 


We may also compare the f spectra of the same partial wave in different mass intervals. In Fig. 35 
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4++1+ p(770) n G 




Figure 34: (color online) The t' dependence of the intensity of the 4++ 1+ p(770) ttG and /2(1270) n F 


waves, integrated over the mass interval around the 04(2040) as indicated by the shaded regions in Figs. 15e 


and 15f The solid red curves represent double-exponential fits using Eq. (461; the fit components are 


shown as dotted curves. See text for details. 
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Figure 35: (color online) The t' dependence of the intensity of the 0 “*■ 0+ [ttttJs tt S wave in two different 
mass regions, which correspond to the two peaks that are indicated by the shaded bands in Fig. |25b[ The 


red curves represent single-exponential fits using Eq. (46) with A 2 = 0. The fitted t' ranges are indicated 


by the solid curves; the extrapolations are shown as dashed curves. See text for details. 


this is shown for the 0 ''O'*' [TrvrJsvrS' wave using the peak regions 1.1 < < 1.3GeV/c^ and 

1.7 < m 3 ,r < 1.9GeV/c^. These mass intervals contain the low-mass part of a potential 7r(1300) 
contribution and the peak region of the 7r(1800), respectively. 

From the above figures, we can see that at low t', all waves have a large single-exponential 
component. In addition, we can distinguish three characteristically different patterns in the t' 
spectra: i) for about half the spectra, the single exponential dominates the full measured t' range; 
a) many waves show larger deviations at higher t', suggesting additional components; in) a 
few waves exhibit a minimum in the t' region between 0.3 and 0.6 (GeV/c)^ (see e.g. Figs. 32b 


and 35a I. The position of such minima is far below the diffractive minima observed in elastic irp 
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scattering. 


Our general ansatz for the description of the observed t' spectra is a sum of two terms, each 
containing an exponential function multiplied by an M-dependent terni^ of the form {t')^ with 
M > 0: 




(46) 


Here, the Ai are real-valued parameters. The above formula is not able to describe the behavior 
of waves that show minima in their t' distribution. 


For each partial wave, one or two specihc Svr mass ranges are selected, which cover known 
resonances. However, the t' dependence of the intensity in these mass ranges still reflects the Svr 
system as a whole in a given partial wave, with both resonant and nonresonant contributions. Bin 
migration effects due to the limited t' resolution of the apparatus are not corrected for. However, 
in the analyzed range the t' resolution is better than 0.02(GeV/c)^ (see Section 2.31, which 
renders the observed t' spectra only slightly shallower than the true ones. 


We perform two kinds of fits: single-slope fits, where the parameter A 2 in Eq. (46) is set to zero 


and double-slope fits, where all four parameters are left free. For cases where the t' spectra exhibit 
more than one component, the range of the single-slope hts is limited to lower t' values. Since 


Eq. (46) is not able to describe the dip structures appearing in some t' spectra, those distributions 


are fit only with a single exponential. Details on the fit results are summarized in Tables and 
that show the ranges in t' and the resulting slope parameters, the intensity ratio of the two 
components within the fit range, and the fit quality in terms of y^/ndf. It should be noted that 
is calculated using the integrals of the model function over the respective t' bins. About half 
of the spectra require a description with two slopes. For spectra that can be described by a single 
slope only, the double-slope fit results in a second component having a very small relative weight. 
For these cases, the values are omitted from Table 


Because of the high precision of the data, statistical uncertainties on the extracted slope parameters 
are negligible and therefore the uncertainties are mostly of systematic nature. The values of the 
slope parameters depend, among other things, on the choice of boundaries of the mass interval and 
the fit range in t'. Given the complex interplay between resonant and nonresonant components, 
which can only be disentangled later at the stage of the mass-dependent ht [33], we have not 
attempted to quantify the systematic uncertainties. We therefore quote the slope parameters 
rounded to two-digit precision and do not give the respective uncertainties. In the figures, the fit 
functions are represented by red curves. For the double-exponential hts, the full t' range from 
0.1 to 1.0 (GeV/c)^ is used. In contrast, the single-exponential hts are performed using narrower 
t' ranges, which are chosen individually for each partial wave and mass region (see Table [^. In 
this case, the ht ranges are indicated by solid red curves, while the extrapolations to the full t' 
range are shown as dashed red curves. In every t' bin, the integral of the ht function, which enters 
the function to be minimized, is shown as a red horizontal line, while the blue line represents 
the data, so that their difference directly indicates the ht quality. For the double-exponential hts, 
the two components are shown in addition as dotted curves: blue for the steep component and 
green for the shallow one. In the following, we shall discuss the observed characteristics for each 
sector. 


JPC = 0-+: We study the 0 waves containing /o(980) and [Trvrjg as isobars. The intensity 
spectrum of the 0 ^ 0'*' tt S wave shown in Fig. |25b exhibits two pronounced maxima and 


differs strongly from the one for the corresponding wave with a /o(980) isobar shown in Fig. 25a 


The higher-lying maximum in the [TTTrJsvrS' decay mode corresponds to the 7r(1800) and exhibits 


'’’'Given by the forward limit of the Wigner D-functions (see Ref. |62p. 
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Table 6: Slope parameters bi from a single-exponential fit to the t' spectra in the given t' ranges. The 
listed TO 3 ^ intervals, over which the intensity is integrated, cover the peak regions of the different partial 
waves. 


Partial Wave 

msn Range 

t' Range 

bi 

X^/ndf 


[GeV/c2] 

[(GeV/c)2] [(GeV/c)-2] 


1++ 0+/>(770) 7r5 

[1.10,1.30] 

[0.100,0.326] 

12 

120 

1++ l+/>(770) 7r5 

[1.10,1.30] 

[0.100,0.724] 

16 

6.6 

1++ 0+/2(1270) 7rP 

[1.68,1.88] 

[0.100,1.000] 

8.4 

4.3 

2++ 1+p{770)7tD 

[1.20,1.40] 

[0.100,0.326] 

8.9 

37 

2++ 2+ p{770)ttD 

[1.20,1.40] 

[0.164,0.724] 

8.5 

18 

2++ l+/2(1270) 7rP 

[1.20,1.40] 

[0.127,0.724] 

7.5 

5.3 

2-+0+p(770) 7rF 

[1.56,1.76] 

[0.113,0.724] 

9.2 

5.1 

2“+ 0+/2(1270) 7r5 

[1.56,1.76] 

[0.100,0.326] 

9.8 

30 

2"+l+/2(1270) 7r5 

[1.56,1.76] 

[0.113,0.724] 

6.3 

3.5 

2"+ 0+/2(1270) 7rZl 

[1.56,1.76] 

[0.113,1.000] 

7.8 

2.7 

2-+0+/>(770) 7rF 

[1.80,2.00] 

[0.113,0.724] 

7.2 

3.5 

2“+ 0+/2(1270) 7rZl 

[1.80,2.00] 

[0.113,0.724] 

8.4 

14 

4++ l+p(770) 7rG 

[1.86,2.06] 

[0.164,0.724] 

8.8 

25 

4++ l+/2(1270) 7rF 

[1.86,2.06] 

[0.164,0.724] 

8.4 

11 


Waves with /o(980) isobar 



0“+0+/o(980)7r5 

[1.70,1.90] 

[0.100,0.724] 

11 

5.6 

1++ 0+/o(980)7rP 

[1.38,1.58] 

[0.100,0.724] 

11 

2.1 

2"+0+/o(980)7rF) 

[1.56,1.76] 

[0.100,0.724] 

8.4 

4.2 

2"+0+/o(980)7rF) 

[1.80,2.00] 

[0.100,0.724] 

7.3 

4.9 


Waves with [Trvrjs isobar 



0 ^ O"'' tt S 

[1.10,1.30] 

[0.100,0.326] 

22 

55 

0 ^ O"'' [7r7r]5 TT S 

[1.70,1.90] 

[0.100,0.449] 

12 

4.1 

O"*" [7r7r]5 TT P 

[1.10,1.30] 

[0.100,0.449] 

13 

11 

2 ^ O"*" [vrvrj^ tt D 

[1.56,1.76] 

[0.100,0.724] 

11 

35 

2 ^ O"*" [vrvrj^ tt D 

[1.80,2.00] 

[0.100,0.724] 

11 

5.8 
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Table 7: Same as in Table but for the double-exponential fit over the full t' range of 0.1 < t' 

< 1.0 (GeV/c)^. The given intensity ratio is the ratio of the integrals I 1 /I 2 of the exponentials with slopes 
bi 2 - Partial waves with t' spectra, which exhibit a clear dip structure (marked with f) or which are already 
well described by a single slope as shown in Table (marked by t) are not fit with the double-exponential 
model. 

Partial Wave 

Range 

[GeV/c2] 

bi 

[(GeV/c)-2] 

&2 

[(GeV/c)-2 

Intensity 
] ratio I 1 /I 2 

X^/ndf 

1++ 0+p(770) 7r5 

[1.10,1.30] 

15 

8.0 

1.2 

6.9 

1++ l+p(770) 7r5 ^ 

[1.10,1.30] 

— 

— 

— 

— 

0+/2(1270) 7rP ^ 

[1.68,1.88] 

— 

— 

— 

— 

2++ 1+ p{770) IT D 

[1.20,1.40] 

12 

6.0 

0.69 

7.3 

2++2+ p{770) IT D 

[1.20,1.40] 

12 

5.7 

1.2 

2.2 

2++l+/2(1270) 7rP 

[1.20,1.40] 

13 

6.5 

0.27 

3.1 

2-+0+p(770) 7rF ^ 

[1.56,1.76] 

— 

— 

— 

— 

2“+ 0+/2(1270) 7r5 

[1.56,1.76] 

14 

7.0 

0.71 

4.1 

2“+l+/2(1270) 7r5 t 

[1.56,1.76] 

— 

— 

— 

— 

2“+ 0+/2(1270) 7rZl t 

[1.56,1.76] 

— 

— 

— 

— 

2-+0+p(770) 7rF t 

[1.80,2.00] 

— 

— 

— 

— 

2"+ 0+/2(1270) 7rZl 

[1.80,2.00] 

11 

5.9 

1.5 

0.47 

4++ l+p(770) 7rG 

[1.86,2.06] 

13 

6.3 

1.1 

4.1 

4++ l+/2(1270) 7rF 

[1.86,2.06] 

14 

6.1 

0.84 

1.9 


Waves with /o(980) isobar 



0"+0+/o(980)7r5' ^ 

[1.70,1.90] 

— 

— 

— 

— 

0+/o(980)7rP t 

[1.38,1.58] 

— 

— 

— 

— 

2-+0+ fo{980) TT D t 

[1.56,1.76] 

— 

— 

— 

— 

2-+0+ fo{980) TT D t 

[1.80,2.00] 

— 

— 

— 

— 


Waves with [ttttJs 

isobar 



0 ^ 0^ [Trvr]^ TT S ^ 

[1.10,1.30] 

— 

— 

— 

— 

0 ^ 0^ [Trvr]^ TT S ^ 
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Figure 36: (color online) The t' dependences of the intensities of the 1++O’*' waves with [ttttJs (a) and 
/o(980) (b) isobars taken around the peak region in the respective wave as indicated by the shaded bands 
in Figs. |25^ and 25f The red curves represent single-exponential fits using Eq. (46) with A 2 = 0. The 
fitted t' ranges are indicated by the solid curves; the extrapolations are shown as dashed curves. See text 
for details. 


a slope parameter of 6 ~ 12 (GeV/c) ^ (see Fig. 35b|, similar to that for the /o(980) tt S decay 


mode. This is in agreement with the expectation that a resonance should have the same slope 
parameter independent of its decay mode. In both cases, the t' spectra are purely exponential. 
In contrast, the t' spectrum corresponding to the broad structure in the 1.1 to 1.3GeV/c^ mass 
range around the elusive 7r(1300) exhibits a pronounced intensity minimum around 0.35 (GeV/c)^ 
and a second maximum around 0.6 (GeV/c)^ (see Fig. 35a). This behavior suggests that different 
production processes are interfering and is similar to predictions by the three-component Deck 
model |61) . The single-exponential fit to the low-t' region results in an exceptionally steep slope 
of 6 PS 22 (GeV/c)“^. The strikingly different t' dependences of the 7r(1300) and 7r(1800) mass 


regions are further illustrated by Fig. 24 and in Section 6.3 


The mass region around the ai(1260) peak contains both resonant and non¬ 
resonant contributions, the latter ones dominated by the Deck process. Using a single slope, 
we obtain h ~ 12(GeV/c)“^ for the ai(1260) region in the I"’"*'0"*“ p(770) vr S' wave, a similar 
value of 6 ~ 13(GeV/c)“^ in the I"*"*'0^ [7r7r]5 7rP wave (see Fig. 36a), and a steeper slope of 
b ^ 16(GeV/c)“^ for the I"'"'' I"''/o(770) vr S wave (Fig. 32b). The t' distribution of the M = 0 
wave is much better described by two slopes with nearly equal intensity (Fig. 32a). The one 
for the M = 1 wave exhibits a dip at approximately 0.5(GeV/c)^. However, the slope bi of 
the steep component in the double-exponential fit of the M = 0 wave is similar to that of 
the M = 1 wave extracted using the single-exponential model in the region of lower f. If we 
interpret the components with the steep slopes to be of nonresonant origin, we would conclude 
that it contributes about 50 % to the M = 0 intensity and that it dominates the M = 1 wave. 
The intensity around the peak in the 1+"*“ 0"*“ /2(1270) n P wave at 1.8GeV/c^ exhibits a nearly 
single-exponential t' spectrum with a slope similar to that of the 02(1320) and 04(2040) regions 
in the respective 2+'*' and waves. The Oi(1420) peak in the I"*"*' 0"^ /o(980) tt P wave is well 
described by a single exponential (see Fig. 36b) and has a slope parameter of 6 ~ 11 (GeV/c)“^ 
similar to that of the 7r(1800) in the 0 waves with /o(980) and isobar. This finding is 

consistent with a slope parameter of 5 ~ 10 (GeV/c)“^ that was extracted for the ai(1420) in a 
mass-dependent fit |55| . 
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= 2++: The waves containing the 02(1320) are best described using two exponentials and 
show similar behavior regardless of the type of the isobar and the orbital angular momentum L in 
the decay. A shallow component with 62 ~ 6 (GeV/c)“^ is accompanied by a steeper component 
of comparable magnitude with bi ~ 12 (GeV/c)“^. Different spin projections M are equally well 
described (see Fig. 33). In the waves, the steep components cannot be directly identified 
with nonresonant contributions, because they are small. It cannot be excluded that the two 
components are caused by the interference of the low-mass tails of excited 02 states with the 
ground state, which may contribute to the t' spectra with different slopes. 


= 2 This is studied in four partial waves containing /9(770) and /2(1270) isobars 
and two waves containing [vrvr]^ and /o(980) isobars. The latter two show striking interference 
effects and are discussed further in Section The t' spectra are studied in two different mass 
intervals: one containing the 7r2(1670), the other the 7r2(1880). The observed pattern is rather 
irregular. Single-exponential fits yield slope values from about 6 to 11 (GeV/c)“^. For the mass 
region around 7r2(1670), the t' spectra fall into two classes: i) distributions that are single¬ 
exponential or have only small contributions from a second slope [2 O'*" p(770) vr F, 2 I"*" 

72(1270) vr S, 2 ^ 0^ 72(1270) vr D, and 2 O'*" 7o(980) vr D waves] and ii) distributions that need 

two exponentials [2 ^ O'*" 72(1270) vr 5 and 2 O"*" [vrvrjgvrD waves]. The latter waves have a 

shallower component with a slope of around 7 (GeV/c)“^ accompanied by a steep component of 
similar magnitude. The slopes of the single-exponential spectra vary considerably. 

The pattern is different for the vr2(1880) mass region. Here, the 2 ^ p{770) tt F, 2 O'*" 

7o(980)vrL), and 2 O"*" [vrvr] 5 vrZl waves are nearly single-exponential. However, the latter has 

a steeper slope of 11 (GeV/c)“^ compared to 61 ~ 7 (GeV/c)“^ for the former two. The 2 
0^ 72 ( 1270 ) TT D wave requires two slopes, where the steep slope is about 11 (GeV/c)“^ and the 
shallow one approximately 5.9 (GeV/c)“^. 

The waves containing the 04 ( 2040 ) are studied in decays into two different isobars. 
The t' spectra follow the pattern observed for with one slope of 62 ~ 6 (GeV/c)“^ 

and a steeper component described by 61 ~ 13 to 14 (GeV/c)“^ of about equal strength. 

In summary, for single-exponential fits of the t' spectra, we find a general trend of shallower slopes 
with increasing mass. Waves with dominant resonant contributions, like e.g. the 2~^^ and 
waves, have slopes in the range from 7 to 11 (GeV/c)“^, which are equal for different decay modes. 
In contrast, waves with large nonresonant contributions, like e.g. the 1"’“"’' p{770) vr waves, show 
typically steeper slopes in the range of 12 to 16 (GeV/c)“^. Many waves are better described 
by a two-exponential model. However, in general the two components do not seem to separate 
nonresonant from resonant contributions. This may be due to possible large interference effects 
or contributions from excited states. Signs of such interferences are observed in the t' spectra of 
some waves, which exhibit a dip around t' ~ 0.3 to 0.6 (GeV/c)^ and thus can be described by 
the single-exponential function only in a limited t' interval. 

Our results of the fits using single exponentials can to some extent be compared to earlier analyses 
done on this topic. The single-slope parameters in the mass region around the 02(1320) of 
7.5 to 8.9 (GeV/c)“^ agree well with the results obtained for the ? 7 vr“ and t]' 7 t~ channels studied 
at the same incident energy. In the former channel, which is dominated by the 02(1320), the 
slope parameter is 8.45 (GeV/c)“^ [50j . In the rj' 7r~ channel, a slope parameter of 8.2 (GeV/c)“^ 
is found in the 02(1320) mass region |63] . As in the present case, natural-parity transfer (M = 1) 
is strongly dominant. Hence, all 02(1320) production characteristics are consistent with being 
independent on the decay channels, as required for true resonances. For O4(2040) production, the 
3vr and r/'vr” results (see Table and Ref. [63], respectively) are consistent with this requirement 
as well. 
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The ACCMOR collaboration |14l [T5l [131CS] pioneered such t' fits for selected waves in the Stt 
mass region between 0.8 and 1.9GeV/c^, describing the t' spectra in the range t' < 1.0 (GeV/c)^. 
For the I"’"*' O"*" p{770) n S wave, which contains the ai(1260), the authors quote an overall slope 
parameter of 6 = (10.1 ± 0.9) (GeV/c)“^, which is similar to our data. For the 2^^ f*" p{770) 
TT D wave, they observe b = (7.3 ± 0.1) (GeV/c)“^ as compared to 6 ~ 8.9 (GeV/c)“^ quoted in 
Table Finally, for the waves 2 0^ /2(1270) vr S and 2 ^ 0"*“ [7r7r]5 vr S, they have extracted 

values of 6 = (8.5 ± 0.3) (GeV/c)“^ and b = (10.7 ± 1.1) (GeV/c)“^, respectively, while selecting 
a mass window from 1.6 to 1.7GeV/c^. Both values are in good agreement with our findings (see 
Table [^. The authors concluded that owing to strong nonresonant effects, the true values for b 
in direct resonance production might be around b = 7.5 (GeV/c)“^. This value agrees with our 
findings of b being in the range of 7 to 11 (GeV/c)“^ for the 2^"*“, d"’"*', as well as for the 0 0^ 
/o(980) TT S and l^^ 0+ /o(980) vr P waves, which we also ascribe to resonant production. 

Results from BNL E852 originate from a mass-independent fit in 12 bins of t' and are shown in 
Ref. |19| for the waves 2 ^ 0^ /2(1270) vr S', 2'’'+ I"*" p(770) vr D, I"*" p(770) vr G, and 1 ^ 1+ 
p{770) vrP for 0.1 <t'< 0.5 (GeV/c)^, but without discussing a functional description of the t' 
dependence in detail. 


6 Determination of tttt S-Wave Amplitudes 


As shown in Section 4.4, the = O"*"^ isobars decaying into vr vr"*" in an S-wave are important 


intermediate states in 3vr meson decays. In the considered 3vr mass range, they consist of i) a 
broad continuum, which is usually described by a parametrization extracted from vrvr S-wave 
elastic-scattering data, and ii) at least two distinct resonances, /o(980) and /o(1500). The much 
debated /o(1370) was not included as a separate isobar in the analysis described in Section]^ 
The key issue is: to what extent the information extracted from vrvr elastic scattering can be 
used to describe spectral shapes and phases of the two-pion O"’"*' isobars in many-body decay 
amplitudes? 

As discussed above, the Oi(1420) appears only in the /o(980) vr P-wave and its strength (but not 
its shape) reveals some dependence on the detailed parametrization used for the /o(980), i.e. a 


Breit-Wigner or a Flatte amplitude (see Section [4.6| and Appendix B.3). This section addresses 
in particular the question whether the observed ai(1420) is truly related to the narrow /o(980) 
or whether it is an artifact of the isobar parametrizations employed in the fit. This is relevant for 
the significance of the new observation as well as for the interpretation of the ai(1420). 


6.1 Method of extracting Isobar Amplitudes from Data 

The conventional isobar model uses fixed amplitudes for the description of the vr^vr'*' intermediate 


states ^ (see Section 4.1). However, we cannot exclude that the fit results are biased by the 


isobar parametrizations used. This is particularly true for the 0"’“"’' vr“vr^ isobars, where we have 
separated a broad vrvr S-wave component from the /o(980) and /o(1500) resonances. In order to 
solve this problem, a novel method inspired by Ref. |64) was developed. It allows us to determine 
the overall amplitude of the O^'*' vr^vr'*' isobars directly from the data. 

For selected isobars, the new method abandons the fixed description of the mass-dependent 
amplitudes which appear in the two-body isobar decay amplitude of Eq. ^ and are part 

of the full decay amplitude of the state X~ defined in Eqs. (11) and (21). The latter amplitude 
factorizes into a part, K,%, which depends on the spherical angles I7 gj = (cosr^GJj <('ty) in tbe 
Gottfried-Jackson frame as well as 17 hf = (cos'duFj </>hf) in the helicity frame, and a second part, 
Aa, that is the mass-dependent isobar amplitude. Taking into account the Bose symmetrization 
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according to Eq. (12), we write for a particular bin 

r23) = AaM- 


(47) 


The two terms represent the two possible tt vr'*' combinations of the vr^ 7^2 vr^ system. The index a 


defined in Eq. (16) represents the quantum numbers of the dvr partial wave. This includes the 


quantum numbers of the tt 7 r+ subsystem 

In our new freed-isobar method, we replace the fixed parametrizations for Aa{m^) by a set of 
piecewise constant functions that fully cover the allowed mass range for m^. The isobar line 
shape is rewritten as: 

^ ^a,k (48) 

k 

where the index k runs over Tr^Tr"*" mass bins. These bins are defined by sets of window functions 
{Uk^^im^)} that are non-zero only in a narrow interval in the isobar spectrum given by the 
bin borders 

1 if <m^< mk+i,^, 


nk,^{m^) = 


0 otherwise. 


(49) 


Here, the bin width 5m^ = — ruk^^ may depend on the mass region of the tt Tr"*" system 

considered. The ^a,k are unknown complex numbers that determine the binned amplitude 
Aa{m^). 


The intensity distribution in a given Svr mass bin, as defined in Eq. (24), contains terms of the 
form 


'7;"l^a(n3,r23) 

k k 

Absorbing the unknown isobar amplitudes £^a,k hrto the transition amplitude via 


T^a.k — ki ^a,ki 


(50) 


(51) 


the bins appear in the intensity via the index k that is summed over coherently, in the same 
way as the partial-wave index a is: 


£=±1 


^(^) = E EE'^A^a,fc(n3,r23) 


where 


Kkiris, t23) = nfi) 77fc,^(mi3) + /c^(i^2^^, nk,^^ 


T 2iflat) 


GJ r)HF^ 


(52) 


(53) 


This means that each 27r mass bin can be treated like an independent partial wave. In this way, 
the same procedure as for the standard mass-independent ht can be used. The hts in Stt mass bins 


yield transition amplitudes 7^ that now depend on m^Tr and 171 ^^-^+. According to Eq. (51), these 


amplitudes contain information on the Svr system as well as on the vr vr^ subsystem. It should 
be noted that the method is restricted to rank 1. Therefore, the rank index was omitted in the 


above formulas. It was discussed in Section 4.2 that rank 1 is sufficient for the positive-reflectivity 


waves. 


In the ansatz for the decay amplitude in Eq. (47), the isobar mass-dependent amplitude Aa{m^) 
depends on the Svr partial-wave index a, i.e. the model permits different isobar amplitudes for 
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Table 8: Borders of the four non-equidistant t' bins, in which the partial-wave analysis with freed isobars 
is performed. The intervals are chosen such that each bin contains approximately 11.5 x 10® events. 

Bin 12 3 4 

t' [(GeV/c)2] 0.100 0.141 0.194 0.326 1.000 


different intermediate states X . This is in contrast to the conventional approach, which uses 
the same isobar parametrization in different partial waves. 


vr' 


The reduced model dependence of the new method and the additional information about the vr” 
subsystem lead to a considerable increase in the number of fit parameters in the mass-independent 
fit. Thus, even for large data sets, the freed-isobar approach can only be applied to a subset of 
partial waves. In the analysis presented here, we replace the fixed parametrizations of the set of 
= O'*”*" isobars, which consists of [TTTrjg, /o(980), and /o(1500), by a set of single piecewise 
constant functions representing the overall dynamical amplitude of all O'*"*" isobars as defined in 


Eqs. (48) and (49). In the following, we shall denote the freed O'*"*' isobar amplitude by [7r7r]o++. 


We determine the [7r7r]g++ amplitudes simultaneously for the waves 0 0^ [7r7r]o+-i- tt S, 0^ 

[7r7r]Q++ ttP, and 2 ^ O'*' [7r7r]Q++ tt D, which are the dominant waves with 0 '’'+ isobars. These 


partial waves replace a set of seven waves with conventional isobar parametrizations (see Table A.l 
in Appendix]^. For all other amplitudes, we keep the isobar parametrizations as discussed in 
Section 4.1 The fits are performed in msTr bins with 40MeV/c^ width, i.e. twice as wide as used 


in the conventional analysis. Each ht results in an Argand diagram for the [7r7r]g++ amplitude 
ranging in the two-pion mass from 2m7r to — ‘m-iT- The bin width in the vr^Tr"’' subsystem 
mass is 40MeV/c^, except for the region 920 < m^-^+ < 1080MeV/c^ around the /o(980)fi^ 
Here, finer bins of lOMeV/c^ width are chosen in order to better resolve the resonance structure. 
In total, 62 two-pion mass bins are used. In order to obtain reasonable statistical accuracy, we 
perform this analysis in only four bins of t', which are listed in Table 

As in the conventional analysis, multiple fit attempts are performed with randomly chosen starting 


values for the decay amplitudes in Eq. (52). Here, the ht with the highest likelihood is 
selected from a set of 50 attempts. For 37r masses below about 1 GeV/c^, the hts tend to be 
trapped in local maxima that deviate from each other only by a few units of log-likelihood. Such 


a behavior is also observed in the hxed-isobar hts (see Section 4.6) 


6.2 Comparison with the Fixed-Isobar Method 

In order to compare the new freed-isobar method with the conventional analysis scheme, the 
hxed-isobar ht was repeated with the coarse binning in 37r mass and t'. Based on this ht, the 
amplitudes of partial waves with the same X~ quantum numbers but different 0"’"’“ isobars, i.e. 
[vrvrjs, /o(980), and /o(1500), are added coherently. For = 0 and 2 ^ of the 37r 


system, the resulting intensities are shown in Fig. 37 as blue data points in two t' regions chosen 
as examples. These spectra are related to those in Fig. which show the intensity distributions 
separately for the [vrvrjs and /o(980) isobars, integrated over the full range of t'. The striking 
interference effects observed in the 2 O"*" /o(980) tt and ['7r7r]5 vr isobaric waves are washed out 

in the coherent sum of the two. 

In our novel approach, we do not separate the different O"*"*' isobar components but obtain the 
overall 0"’"’“ amplitude in bins of and m^-K- This also implies that the correlation of the 

relative phases between the components across the spectrum is not predetermined. The 

red data points in Fig. 37 show the 37r mass spectra for = 0 1"’"'', and 2 obtained 


'‘^Also the first m^-^+ bin, which covers the mass range from 2m.,r to 320MeV/c^, has a slightly different width. 
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Figure 37: (color online) Intensities of the coherent sum of partial waves with the same quantum numbers 
but different 0+“'' isobars, as obtained in the conventional analysis scheme with fixed isobar amplitudes 
(blue), and the corresponding intensities obtained from the freed-isobar fit (red). Left column: low t'; 
Right column: high t'. 
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by coherently summing over all two-pion mass slices [represented by index k in Eq. (48)] in the 


two chosen t' bins. These intensity distributions can be compared directly to those obtained by 
coherently summing over the O'*"*' isobars using the conventional analysis method shown as blue 


data points in Fig. 37 The agreement is good in general, in particular the 7r(1800) region in the 
0 wave matches well. In the wave, the region < 1.2GeV/c^ is enhanced in the fit 
result for the freed isobars except for the highest t' bin. This is partly due to the fact that waves 
with freed isobars have no Svr mass thresholds in the new fit, whereas in the conventional fit. 


m 


the O"'' /o(980) vrP wave was limited to the region of > 1180MeV/c^ (see Table A.l 
Appendix | a|. The largest differences appear in the 2 ^ wave, where we observe in the freed-isobar 
fit an increased intensity in the region around the 7r2(1670) across all t' bins. Systematic studies 
indicate that imperfections in the description of the other isobars used in the PWA ht have an 
influence on the O"*"*' sector. The agreement between the results of the two methods validates the 
parametrizations of the tttt S-wave isobars that are employed in the simpler fixed-isobar fit. 


6.3 Correlation of 27r and Stt Mass Spectra for freed tttt S'-Wave Isobars 

It is instructive to look at the correlation of the [7r7r]Q++ mass spectrum with the Svr mass spectrum 


39 and 1401 are extracted from the 


in different partial waves. The examples shown in Figs. 38 
0 O'*" [7r7r]o++ vr S, 1"’"’“ 0"*“ [7r7r]o++ tt P, and 2 0"^ [7r7r]o++ tt D waves, respectively. The z axis 

of the two-dimensional representations (left columns) is given by |7^^(m^-^+, m37r)P, which is 
normalized such that it can be interpreted as the number of events per unit in 771^-,^+. The 
apparent dependence of the shape of the 2 tt mass distribution on 1713^^ and on of the 37r 
system reveals the different coupling of 37r resonances to the various components of the 27r 
subsystem. In the following, we will discuss the features for each three-pion in detail. 

JPC = 0-+: In the conventional analysis, the coherent sum of the waves with fixed O"’"*' isobars 
exhibits two peaks in the intensity that may be identified with the 7r(1300) and the 7r(1800). 


These two peaks appear very similar in the fit with freed O"’”’' isobars (see Figs. 37a and 37b). 
Since all three pions are in a relative S-wave in the 0 0"^ [7r7r]o++ vr S wave, it is very sensitive 

to nonresonant 37r contributions. In particular, in the wave with the fixed [vrvrj^ isobar, the region 
around = 1.3GeV/c^ seems to have nonresonant components (see Figs. 24, 25b and 35a). 


Also in the freed-isobar fit, the shape of the [7r7r]o++ intensity and its considerable t' dependence 


in the 7r(1300) region suggest this to be mostly nonresonant (see left column of Fig. 38). These 
observations are in accordance with quark-model calculations for the first radial excitation of the 
pion |51) which predict a strong suppression of the [vrvr]^ tt decay mode as compared to /?(770) vr. 
Apparent enhancements of this wave in the 7r(1300) region, which were observed in diffractive pion 
scattering by the VES and BNL E852 experiments [201 Ei, are consistent with our observations 
and were attributed to the Deck mechanism mm- 

In order to study the role of /o(980) for the 0 wave, we sum up the intensity in the mass region 
0.96 < mjr-7r+ < l.OOGeV/c^, which contains almost the full /o(980) and which is indicated by a 
pair of red horizontal lines in the left column of Fig. 38 The resulting 1713 ^ intensity spectra are 


shown in the central column of Fig. 38 for the four t' bins. This simple method does not take into 
account the interference of the /o(980) with the broad Trvr S-wave component. The separation of 
amplitudes would only be possible by fitting the m^-^+ and tustt dependences of the amplitudes. 
The msTT intensity distributions exhibit a clear signal for the 7r(1800). In contrast, no clear 
correlation with a possible 7r(1300) can be identified. We observe the low-mass structures in msjr 
to vanish with growing t'. This indicates the existence of considerable nonresonant contributions 
in this 37r mass region. In a similar way, the role of the /o(1500) is investigated by summing the 
intensity over the range 1.44 < 771^-,^+ < 1.56GeV/c^ as indicated by a pair of blue horizontal 
lines in the left column and shown in the right column of Fig. Again, a clear correlation with 
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Figure 38: (color online) Intensity of the [ 7 r 7 r]Q++ tt 5-wave component of the = 0 ^ 0+ amplitude 

resulting from the freed-isobar fits in four t' bins. Left column: two-dimensional representation of the 
intensity of the 0 ^ O’*" [7r7r]Q++ tt S wave as a function of m^-^+ and Central and right columns: 

intensity as a function of 7723 ^ summed over selected m^-^+ intervals around the /o(980) (center) and 
around the /o(1500) (right) as indicated by pairs of horizontal dashed lines in the left column. The vertical 
dashed lines indicate the centers of the 772 . 3 ^ bins discussed in Section |6.4| 
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Figure 39: (color online) Same as Fig. 
column shows the m^-^+ interval around the /o(980). 


but for the l+“'' O’*" [7r7r]o++ tt P wave, except that the right 
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the 7r(1800) is observed. 


= 1++: The intensity correlations shown in the left column of Fig. 39 are dominated by 
a broad maximum between 0.6 and 0.8GeV/c^ in For increasing t', this structure shifts 

from ttistt = 1.2 to 1.4GeV/c^, almost reaching the ai(1420) region for the highest t' bin. This 
behavior suggests the existence of large nonresonant contributions, which obstruct the observation 
of a possible coupling of the ai(1260) to the broad component of the vrTT S'-wave. The right 


column of Fig. 39 shows the intensity summed over the /o(980) mass region of 0.96 < m^-^+ 
< l.OOGeV/c^. For all t' bins, the mass spectra show a clear ai(1420) peak and no contribution 
of the ai(1260). This demonstrates that the observed ai(1420) signal in the /o(980) vr channel is 
not an artifact of the 0 ^"*“ isobar parametrizations used in the conventional analysis method. 


= 2“^': The 2 ^ intensity correlations shown in the left column of Fig. 40 exhibit a 
vertical band around = 1.5GeV/c^ [below the 7r2(1670)]. The m^-^+ distribution peaks 
below the /o(980). This structure changes its shape and relative strength with t'. The role 
of the /o(980) and /o(1500) isobars is again illustrated by summing the intensities over the 
respective 27r mass ranges, which are shown in the central and right columns of Fig. 40 For 
both 7 r“ 7 r^ mass regions, we observe a clear signal for the 7r2(1880). The intensity maximum 
around 7713 ^ = 1.6GeV/c^ in the /o(980) slice changes its shape and position with t', and hence 
looks different from the 772 ( 1670 ) peak, as it is, for example, observed in the /2(1270) tt decay 
mode and shown in Figs. |18a 18c 


22a, and 22c In addition, the position of the peak is t' 


dependent, which indicates a nonresonant contribution. The VES experiment |20j has reported on 
an excited 1^2 resonances at 2.09GeV/c^. It was observed as a 520MeV/c^ broad enhancement 
in the [ttttJs tt and /o(980) vr waves. We also observe a similarly broad structure in the 37 r system 
at m 37 r ~ 2.2GeV/c^ correlated with a broad bump at 277 masses of approximately 750MeV/c^. 
The shape in seems to change as a function of t'. With the present analysis, we cannot 

confirm the resonance interpretation of this structure. While a corresponding peak is observed 
around m 37 r = 2.2 GeV/c^ in the conventional fixed-isobar fit in the 2 O'*" /o(980) tt D wave (see 


Fig. 25c), no pronounced correlation with the 277 system in the /o(980) mass region is seen in the 
freed-isobar fit. 


6.4 Argand Diagrams and 27r Mass Spectra for freed tttt S-Wave Isobars 

The previous section discussed mainly the correlation of the tttt S-wave and the 377 partial-wave 
intensities. The two-dimensional transition amplitudes extracted from the data furthermore 
contain information about the 7713 ^ and the m^-^+ dependences of the relative phases. These 
phases are measured with respect to the 1"’“'*' 0^ p(770) tt S anchor wave as a function of 277 
mass. They give insight into the composition of the 7777 S-wave amplitude. In order to study 
the influence of the 377 system, we look at the 277 invariant mass spectra for three 377 mass bins, 
chosen below, at, and above clear 377 resonance signals. The centers of the 7713 ,^ bins are indicated 


by green vertical lines in the left columns of Figs, p^to 40 


tPC 


0 The wave with the freed [7777)0++ isobar shows a clear signal for the 7r(1800) 


coupling to /o(980) 7 r and /o(1500) 7 r. The left column of Fig. 41 shows the [ 7777 ) 0 ++ intensity 
as a function of m^-^+ for three different values of 7773 ^, i.e. below, at, and above the 7r(1800) 
resonance for the region of larger t', where the resonance signal is clearer. The three m^-^+ 
spectra are similar. We observe prominent signals for /o(980). Because of phase space, the 
/o(1500) peak appears only in the two higher 377 mass intervals. The enhancement of both states 
relative to the broad tttt S-wave component is strongest at the 7r(1800) mass. The right column 
of Fig. 41 shows the corresponding Argand diagrams measured with respect to the 0"*“ p{770) 


77 S wave. For a fixed mass of the 377 system, the Argand diagram describes magnitude and phase 
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of the 27r amplitude. The phase of the [ 7 r 7 r]g++ amplitude spans nearly two full circles about 
the origin in the 27r mass range from threshold to about 1.6GeV/c^. This reflects the resonance 
character of the /o(980) and /o(1500). The positions of /o(980) and /o(1500) (marked by the blue 
line segments in the Argand diagrams) rotate counterclockwise with increasing 77137 ^, reflecting 
the growing phase of the 7r(1800) with respect to the anchor wave. We conclude a clear coupling 
of 7r(1800) to both /g states, which is more pronounced than that to the broad component of the 
TTTT 5-wave. There is no evidence for the /g(1370) in this wave. The observed behavior of the 
[ 7 r 7 r]g++ phases corroborates the conclusions drawn from the intensity correlations in Fig. 38 In 


the 7r(1800) region, the gross features of the [ 7 r 7 r]g++ phase motion are similar to those observed 
by the BABAR experiment in the 27r subsystem of Df —)• 7 r'*' 7 r“ 7 r+ decays (65j . Differences in 
details are probably due to different nonresonant contributions in the two processes. 

= 1"*"“*": In this wave, the ai(1420) region is most interesting. The extracted [ 7 r 7 r]g++ 
intensities are shown in the left column of Fig. |^for three different values of in a similar 
way as explained above. For 7773,7 values around 1.4GeV/c^, a signal for the /o(980) appears 
sitting above a broad tttt 5-wave structure. The complicated shape of the Tr^vr"’' amplitude is 


again illustrated by the Argand diagrams shown in the right column of Fig. 42 Here, the /g(980) 


contribution can be identified by the semicircle-like structure with a shifted origin, which appears 


at the ai(1420) resonance. Comparing Figs. 42d and 42f a significant counterclockwise rotation of 
this structure by about 90° with respect to its center is observed above the ai(1420). This confirms 
the resonance interpretation of the ai(1420) in the /g(980)7r decay. For 777 „-„+ < 0.8GeV/c^, 
the amplitude does not change much with respect to 7773 , 7 . As a consequence, the relative phase 
of the /g(980) component with respect to the broad mr 5-wave structure also changes by about 
90°. Hence the interference pattern changes from partly constructive at the ai(1420) to partly 
destructive at the higher 7773,7 bin. This results in a sharp drop of the intensity in the 27r mass 
spectrum above the /o(980), which is followed by low intensity at higher 777,7-„+. 

= 2 The clearest signal in this wave is the 772 (1880) that couples to /o(980) tt and 
/g(1500) vr. We again study the 27r subsystem in three 37 r mass regions. The intensity distributions 


and Argand diagrams are shown in Fig. 43 Also here we find semicircle-like structures with shifted 
origin that correspond to the /g(980), which is observed as a clear peak in the corresponding 
m^-^+ distributions. At higher values of 111 ^-^^+, we observe an indication of another small 
circular structure in the Argand diagram, which is correlated with a rise of the intensity in the 
277 mass distribution attributable to /g(1500). The shape of the broad vrvr 5-wave component at 
low 277 masses is somewhat different from that of the other waves, exhibiting more intensity close 
to the threshold. 


7 Summary and Conclusions 


We have presented a detailed partial-wave analysis based on the world’s currently largest data 
set of the exclusive 7 r“ 7 r“ 7 r^ hnal state from diffractive pion scattering off a proton target. 
The PWA was performed independently in 100 bins of the 377 mass 7773,7, each of which was 
subdivided into eleven slices of the reduced four-momentum transfer squared t'. We refer to this 
as mass-independent fit. It is based on the largest wave set used so far in a PWA of this hnal 
state, which contains in total 80 waves with positive rehectivity, 7 with negative rehectivity, and 


in Appendix [A|) . In this paper, a subset of 18 partial waves with positive rehectivity accounting 
for 75.8 % of the total intensity was studied in detail. 

From the study of the general characteristics of partial-wave intensities, two classes of waves 
can be identihed: for some waves, the shape of the mass spectrum shows little dependence on t' 


one incoherent isotropic wave representing three uncorrelated pions (see Section 4.2 and Table A.l 
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Figure 41: (color online) The freed [7r7r]o++ amplitude in the 0 ^ O^ [7r7r]Q++ tv S wave in the range 
0.194 < t' < 0.326 (GeV/c)^ for three intervals in below (top row), at (center row), and above 
the 7r(1800) (bottom row). Left column: intensities as a function of m,r- 7 r+i Right column: Argand 
diagrams. The crosses with error bars are the result of the mass-independent fit. The numbers in the 
Argand diagrams show the corresponding mass value of the 7r“7r“'' system. The data points are connected 
by lines in order to indicate the order. The line segments highlighted in blue correspond to the m,r- 7 r+ 
ranges around the /o(980) from 960 to lOOOMeV/c^ and, if phase space permits, around the /o(1500) 
from 1400 to 1560 MeV/c^. The 27r mass is binned in lOMeV/c^ wide intervals around the /o(980) and 
in 40 MeV jc? wide slices elsewhere. The phase of the Argand diagrams is fixed by the !+■*■ O^ p(770) tt S 


wave. 




















66 


xlO® 




W { 717 Z \ ^^71 P 


> 

(U 

O 


0.3 


0 . 2 - 


0.1 


0.194 < f'< 0.326 (GeV/cr 
1.26 < m3^< 1.30 GeV/c2 


+ 


(a) 


0.5 


xlO® 


1.5 2 

[GeV/c2] 

I'^’^O’^ TlP 


> 

(D 

O 


c 

(U 


0.3 


0.2 


0 . 1 - 


0.194 < f'< 0.326 (GeV/cr 
1.38 <m3^< 1.42 GeV/c2 




+ + ,+i 


VII 


(c) 


0.5 1 


1.5 2 

[GeV/c2] 



400 


> 

(U 

o 


200 


B 


-200 

(d) 



0.194 < f 

< 0.326 (GeV/V 


1.38< m 3 ^< 1.42GeV/c2 


0.9% 



1.07^ 


- 

i.lio 

Vo.965 



1^25 

- 


> 


0.299 

_■_■_^_1_^_ 

r_._I_ I_I_^_ 


200 400 600 

ReCr) [(Events/(GeV/c2))‘'^] 


xlO'’ 


> 

u 

O 


c 

(U 


0.3 


0.2 


0.1 


TVP 


0.194 < f'< 0.326 (GeV/V 
1.50 < m3^< 1.54 GeV/c2 


++ 

++ \ 


r** 


(e) 


0.5 


1.5 2 

[GeV/c2] 


300 


-200 

(f) 


RW { 717 Z \ ^^71 P 


0.194 < f'< 0.326 (GeV/V 
1.50< m3^< 1.54 GeV/c2 


0.965 



100 200 300 400 

R&{T) [(Events/(GeV/c2))‘'^] 


Figure 42: (color online) Same as Fig. 
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Figure 43: (color online) Same as Fig. 41 but for the 2+0+ [7r7r]Q++ tt D wave in three bins around 
the 712(1880). 
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(see e.g. Figs. 17b and 17d), while for others moderate (see e.g. Figs. 18a and 18c) or even large 
variations (see e.g. Figs. 17a and 17c I are seen. These variations indicate the presence of weak or 
strong nonresonant contributions that do interfere with the resonant components and may have a 
characteristically different dependence on t'. Among the stable waves, where the peak positions 
do not significantly depend on t' , we find: 


0“+0+/o(980) 7r5 
l++0+/o(980) 7rP 
2++ l+p(770)7rT) 
2++ 2+ p{770)ttD 


2“+l+/2(1270) 7r5 
2“+0+/2(1270)7rT) 
4++ l+p(770)7rG 
4++l+/2(1270) 7rF 


The following waves show significant peak shifts or large variations of the shapes of their mass 
distribution as a function of t': 


0 O"'' [7r7r]5 TT S 

I"’"'*' O'*' [7r7r]5 TT P 
1++0+p(770) TT S’ 
1++1+p(770) vr S 
l++0+/2(1270) 7rP 


2++ l+/2(1270) 7rP 
2 ^ O"'' vr D 
2-+0+p(770)7rF 
2“+0+/o(980) 7rT) 
2“+0+/2(1270)7rS 


In the subset of 18 waves, clear resonance peaks are found in the partial-wave intensities of the 
following decay modes: 


7r(1800) —)• /o(980)7r S-wave 
7r(1800) —>• [7r7r]5 7r S-wave 
ai(1260) —)■ p{770)7r S-wave 
ai(1420) —)• /o(980) tt P-wave 
02 ( 1320 ) —)• p(770)7r D-wme 


7r2(1670) —)■ 72 ( 1270 ) vr S-wave 
7r2(1880) —> 72 ( 1270 ) vr D-wave 
04 ( 2040 ) —)• p{770) TT G-wave 
04 ( 2040 ) —)• 72 ( 1270 ) vr F-wave 


The new ai(1420), which was presented in Ref. [55], is only seen in the l^"*" 0^ 7o(980) tt P wave. 
No evidence for a corresponding resonance structure is observed in I"*"*' 0^ [vrTrj^vrP, nor in waves 
containing other isobars. The 7r2(1670) and 7r2(1880) seem to have different couplings to various 
decay modes. A peak attributable to the 7r2(1670) appears, for example, dominantly in 2 O"*" 

72 ( 1270 ) vrS, while a 7r2(1880) peak is dominant in 2 ^ O"*" 72(1270) tt P. Both states seem to 
couple to the 2 ^ O"*" p{77Q) tt F wave, however, with different apparent strength as a function of 
t' . In turn, only the 7r2(1880) shows a clear coupling to the 7o(980) tt decay mode. The shape of 
the structure observed in this decay mode around ttistt = 1.6GeV/c^ changes as a function of t' 
(see Section 6.3). 


t' Dependences We have investigated the production characteristics of the 7 r“ 7 r“ 7 r'’' final state 
by studying the t' dependence for the overall data sample as a function of 772 , 3 ^ as well as for 
individual partial-wave intensities in 37 r mass regions around known resonances. The fits to the 
overall t' spectra require two exponential functions in order to describe the fall-off with t' (see 
Figs. 30c and 30d). The slopes of the two exponentials and their relative contributions change 
with increasing 37 r mass, leveling off for 7723 ^ > 1.3GeV/c^ (see Fig. 31). 


The slope parameters for individual waves in 37r mass regions around resonances exhibit a complex 
pattern. Qualitatively, mass regions with strong nonresonant contributions are characterized 
by a steep drop-off with t' and thus larger values for the slope parameter up to 22 (GeV/c)“^. 
Gonsiderable deviations from the single-exponential behavior are observed for mass regions 
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around the 7r(1300) in the 0 O"*" [Tr'/rj^vrS' wave and around the Oi(1260) in l+p(770)7r5. 

In these two waves, we find a minimum of the intensity at values of t' of about 0.4(GeV/c)^ 
(Fig. 35a) and 0.6(GeV/c)^ (Fig. 32b|, respectively, which may be attributed to interference 
effects of different production processes. Other distributions can be described well by only a 
single exponential, a parametrization also employed for the ht of the low-t' region in the case of 
dip structures. Mass regions dominated by resonances show typically a shallower drop-off with 
slope parameters between 7 and 11 (GeV/c)“^. However, these regions are often better described 
by a double-exponential model. Hence the observation of a steep component does not exclude a 
dominant resonant contribution. 


For mass regions with clear resonance signals, e.g. 02(1320) and 04(2040), slope parameters 
are found to be similar for different waves belonging to the same even with different spin 

projections M. We have studied the production of waves with = 1"’“'*', 2 and with 
different M. We observe a reduction in their production rate by about an order of magnitude 
with every unit of M. For the 02(1320), the intensity ratio for the two spin projections is in good 
agreement with the one observed in the r]'K~ decay channel |50]. At the same time, we conhrm 
that the t' dependences follow the theoretically expected suppression factor (t')^ at small values 
of t' for M = 0, 1, and 2 (see Figs. 32 and 33). This observation points to the spin characteristics 


of the Pomeron exchange, which is dominant here. 


The l^"*" O'*" /o(980) vrP wave is of particular interest. In the mass region of the new ai(1420), this 
wave exhibits a nearly exponential t' spectrum with a slope parameter of about 11 (GeV/c)“^, 
which is similar to that of the 7r(1800) in the same decay mode. This supports the resonance 
interpretation of the ai(1420) signal. The slope is in agreement with the slope parameter of about 
10 (GeV/c)“^ that was extracted for the ai(1420) in a mass-dependent ht [55]. 

The 2 waves show no clear pattern. Single-exponential hts in different decay modes around 
the '7r2(1670) and 7r2(1880) give slope values between 6 and 11 (GeV/c)“^. For those 2 waves 
that are better described by two exponentials, the dominant slope has a similar range. 


TTTT S-Wave Amplitudes For the hrst time, a detailed study of the amplitude of the Tr^Tr"*" 
S-wave isobar with = O'*' O'*'^ in the decay of the 7r“7r“7r+ system was performed. This 

was achieved by using the freed-isobar technique (see Section]^. The 27r amplitudes are extracted 
independently for different 3'7r partial waves in each bin of and t'. We have presented the 
correlations of the intensities of the independent freed TT~Tr~^ isobar amplitudes with those of 
the 37r system for 0 1"'"'', and 2 ^ three-pion quantum numbers. These correlations 

reveal a selective coupling of 37r resonances to the scalar isobars /o(980) and /o(1500) and less 
clear correlations with a broad Trvr 5-wave component. The new method does not only yield the 
two-dimensional intensity distribution, but also provides information about the full 27r amplitude 
for each msTr bin. In the corresponding Argand diagrams, signals for /o(980) and /o(1500) show 
up as semicircular structures with rapid counter-clockwise motion with increasing In the 

three waves studied, there is no evidence for a distinct /o(1370) resonance in the 7r“7r^ subsystem. 

For = 0 and l^'*' of the 37r system, the mz-w spectra connected to the broad component of 
the TTTT 5-wave show enhancements around msTr = 1.2GeV/c^, which might naively be interpreted 
as ai(1260) and 7r(1300), respectively. These structures significantly change their shape as a 
function of t', thereby suggesting that they are influenced by nonresonant processes. For the 
7r(1800), we observe a coupling to /o(980) 7r and a somewhat weaker one to /o(1500) 7r. The 
Argand diagram shows clear semicircular structures corresponding to the /o(980) and /o(1500). 
Similarly, in the 2 ^ 37r wave, the coupling of the 7r2(1880) to /o(980) tt and /o(1500) vr is seen. 
For the 37r wave with we observe a clear correlation of the /o(980) isobar with the 

new ai(1420) resonance |55| in all bins of t'. This is in contrast to the broad component of the 
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TTTT S'-wave, which shows a strongly t'-dependent correlation with msTr and a shift of the intensity 
maximum towards higher values of mzn with iucreasiug t'. A possible explauatiou of this shift 
is the Deck process. At large values of t', the rapidity gap betweeu the [Trvrj^ system aud the 
bachelor piou is iucreased (see Fig. 19), which leads to higher Stt masses. The shift of iuteusity 
across the Svr mass spectrum with t' could explaiu the complicated behavior of some t' spectra 
(see Sectiou]^. 


Based ou the aualysis described iu this paper, we extracted the properties of resouauces aud of 
uouresouaut coutributious as well as their production characteristics, which will be described iu 
detail iu a forthcoming paper |33) . 
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A Wave Set 


Table [A!T] lists the wave set used for the mass-independent fit. Note that in the reflectivity basis 
= (even)'*'^ waves with = O"*" are mathematically forbidden [see Eq. (19) on page 16 . 


B Systematic Studies of Partial-Wave Analysis Model 
B.l Rank of Spin-Density Matrix 


As pointed out in Section 3.3 ranks > 1 of the spin-density matrix provide a way of modeling 
incoherences between partial waves. This is done by introducing additional sets of transition 
amplitudes. These sets are assumed to correspond to different noninterfering production processes, 
each with its own production phase. By performing the analysis in bins of t', it was found that 
Nr = 1 is sufficient for positive-reflectivity waves. This also leads to higher stability of the 
mass-independent hts. 


In Fig. B.l, we show in red the intensities of selected partial waves obtained from fits with rank 2 
for the positive and negative-reflectivity waves. This is compared to the standard ht (blue data 
points), where rank 2 was used only for waves with e = — 1. In the rank-2 fit, the flat wave 
disappears completely and the intensity of the negative-reflectivity waves is approximately halved. 
Slight modifications of the shape of resonance structures are observed in some partial waves. 
Several partial waves exhibit artificial peak structures in the 1.0 to 1.3GeV/c^ 3'7r mass region, 
like e.g. shown in Fig. B.ld[ Altogether, we prefer to use rank 1 for the positive reflectivity waves. 


B.2 Omission of Waves with Negative Reflectivity 

The PWA model defined in Section has two types of incoherent contributions, rank and 
reflectivity [see Eq. @]- The latter one is determined by the naturality of the exchange particle 
(Regge-trajectory) mediating the scattering process. Including e = —1 partial waves, we effectively 
allow for the exchange of Reggeons other than the Pomeron, which is expected to be suppressed 
at beam energies of 190 GeV. In our PWA model, we have included seven waves with negative 
reflectivity (see Table A.l in Appendix]^. Negative and positive-reflectivity waves have different 
angular distributions. In order to study how well the fit is able to separate the two sectors, we 
have performed fits without any e = —1 waves. The result is shown in red in Fig. B.2 for two 


selected waves. With the exception of the flat wave, the intensities of all waves stay practically 
unaltered. This demonstrates that the positive and negative-reflectivity sectors are well separated 
by the analysis method. 


B.3 Variation of the Isobar Parametrization 


In the employed PWA method, the X~ decay amplitudes ^^9- (24)] are not allowed 

to have any free parameter. Fixed parametrizations for the isobar amplitudes [see 

Eq. ([^] have to be used, which are taken from literature (see Table |^. While eventually these 
parametrizations could be extracted from our data following the analysis scheme outlined in 
Section [^ for this paper we still use the conventional approach. 


The /9(770) is the dominant isobar. As discussed in Section 


4.1, different Breit-Wigner parametriza¬ 


tions exist for the /o(770). Using Eq. (32) instead of Eq. (|® ) for the mass-dependent width F (m) of 
the /9(770) changes the intensity of the structure in the ai(1260) mass region in the I"*”*" O"*" [Trvrjs 
■n P wave and that of the 02(1320) signal in the I"'' /2(1270) tt P wave (see Fig. B.3). Both 


structures seem to be contaminated by model leakage from the respective dominant p{770) tt decay 
modes. The other 16 waves listed in Table [^remain practically unchanged. When using in addition 
the PDG averages for the p{770) parameters of rriQ = 775.26 MeV/c^ and Fq = 149.1 MeV/c^ |9], 
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Table A.l: Wave set used for mass-independent fit: 80 waves with positive reflectivity, 7 with negative, 
plus an incoherent isotropic wave. 


JPC 

Isobar 

L 

Threshold [MeV/c^] 
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Figure B.l: (color online) Comparison of f'-summed partial-wave intensities obtained from the standard 
mass-independent fit with rank Nr = 1 of the spin-density matrix for waves with e = -fl and rank 2 for 
the ones with e = — 1 (blue/black) with the intensities from a fit with rank 2 for both sectors (red/gray). 
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Figure B.2: (color online) Comparison of f'-summed partial-wave intensities obtained from the mass- 
independent fits with (blue/black) and without (red/gray) e = —1 waves. 

the log-likelihood values, summed over the 11 f' bins, decrease by more than 1000 units in the 
mass range between 0.95 and 1.35GeV/c^. 
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Figure B.3: (color online) Comparison of f'-summed partial-wave intensities obtained from mass- 
independent fits using two different parametrizations for the p(770) isobar amplitude: a Breit-Wigner 
with Eq. (40) (blue/black) and one with Eq. (32) (red/gray). 


We have also investigated the sensitivity of the PWA result with respect to changes in the 
parametrization of the /o(980) and ['7r7r]5 isobars. For the /o(980), the Flatte parametrization 
used in the standard analysis is replaced by a modihed S-wave Breit-Wigner amplitude of the 
form: ^ 

^/o(98o)("i;mo,ro) = 312—z:5—(b-l) 


where 


r(m) = ro^. 

Qo 


rriQ — m? — i mo r{m) ’ 
Q 


(B.2) 


The /o(980) parameters are mo = 980MeV/c^ and Tq = 40MeV/c^. Figure B.4 shows in red 


some selected partial-wave intensities from this study. The Breit-Wigner parametrization has 
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less pronounced tails and covers a narrower 27r mass range. This leads to nearly a factor of 
two lower intensities in the /o(980) tt partial waves. The shapes of the resonance structures in 
these waves remain unaltered. Interestingly, the 7r(1800) peak in the 0 O'*' [7r7r]5 7r5 wave also 

decreases when the Breit-Wigner parametrization is used for the /o(980). On the level of the 
mass-independent ht, this behavior cannot be explained. In contrast to the 7r(1800) peak, the 
structure in the 7r(1300) region remains unaltered. The fit with the Flatte parametrization has a 
higher likelihood than the fit with the Breit-Wigner one. 
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Figure B.4: (color online) Comparison of t'-snmmed partial-wave intensities obtained from mass- 
independent fits using two different parametrizations for the /o(980) isobar amplitude: Flatte parametriza¬ 
tion [Eq. ([4^, blue/black] and modified S'-wave Breit-Wigner [Eq. (B.ll with Eq. (B.2|, red/gray]. 


Also for the broad component of the vrvr S-wave amplitude various parametrizations exist. In 
addition to the modified M solution from Ref. sa, we tried the Ki solution from Ref. m with 
the /o(980) pole subtracted, using the modified S-wave Breit-Wigner amplitude of Eq. (B.l) with 
Eq. (B.2) [20]. In order to be consistent, the same /o(980) amplitude was also used for the partial 
waves with the /o(980) isobar. The result is very similar to the one of the fit with the modified 
M solution for the [vrvr] 5 and the S'- wave Breit-Wigner amplitude for the /o(980) isobar discussed 
above. 
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B.4 Variation of Event Selection 


In order to study the potential influence of backgrounds from kaon diffraction, kaon pairs in the 
final state, and central-production reactions, the mass-independent fit was performed on a data 
sample, in which the particle identification in the beam and spectrometer was not used and the 
rejection of central-production events as described in Section 2.3 was not applied. Therefore, 
possible background contributions are expected to be enhanced in this data sample. The effect is 
shown in Fig. B.5 for selected waves. The data sample with looser cuts contains approximately 
20% more events in the analyzed range of and t'. Hence, the partial-wave intensities are 
larger and typically scale proportionally to the sample size. The peak shapes of resonances are 
in general unaffected by the different event selection. The most noteworthy effects of the looser 
cuts are an over-proportional increase of the flat-wave intensity by nearly a factor of two and an 
enhancement of structures at large Svr masses in some waves. 
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Figure B.5: (color online) Comparison of f'-g^mmed partial-wave intensities obtained from mass- 
independent fits for standard event selection (blue/black) and a looser event selection (red/gray), where 
particle identification and central-production rejection have not been applied. 
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B.5 Variation of t' Binning 


In order to study, whether the chosen binning in the four-momentum transfer squared t' has any 
effect on the partial-wave analysis, the t' bins defined in Table and shown in Fig. 11 where 
halved, yielding in total 22 bins. The finer t' binning has practically no effect on the partial-wave 


intensities (see Fig. B.6). Only the flat wave has lower intensity over the full msjr range in the 
case of hirer t' bins. 
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Figure B.6: (color online) Comparison of t'-summed partial-wave intensities obtained from mass- 
independent fits using two different t' binnings: 11 t' bins (blue/black) and 22 t' bins (red/gray). 


C Acceptance 


In the following, we describe the 7r“7r“7r'*' detection efficiency of the COMPASS apparatus in 
absolute terms. For this we have generated Svr events distributed isotropically in phase space 
and passed them through the COMPASS detector simulation and reconstruction chain. The 
same selection cuts were applied as used for the real data. For hxed values of msir and t', the 
acceptance is a hve-dimensional function, of which we show only projections. 


Figure C.l shows the 7r“7r“7r''“ detection efficiency as a function of the two angles, cosi?gj and 
4>ty, of the isobar in the Gottfried-Jackson frame (see Section 3.1 for the dehnition). The 
acceptance is shown in four regions of and t'. The Monte Carlo data show a rather flat 
acceptance with a small dip for in-plane events at forward angles. This structure exhibits some 
dependence on t' and mass. The corresponding distributions in the helicity frame are shown in 
Fig. C.2 Here, no signihcant structures are visible. 
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Figure C.l: (color online) Detection efficiency for 7r“7r“7r+ phase-space events in different regions of 
m 37 r and t'. Each graph shows the detection efficiency as a function of the angles, cosdcj (abscissa, from 
— 1 to -1-1) and 0 ty (ordinate, from —180° to -1-180°), of the isobar in the Gottfried-Jackson frame. 























































t' [(GeV/c)^] 


79 



Figure C.2: (color online) Same as Fig. C.l but for the detection efficiency as a function of the angles, 
cosz^hf and (/>hF) of the 7r“ in the helicity frame. 
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D Additional Partial-Wave Intensities 


In this section, we present the intensities of the remaining 69 waves of the 88-wave PWA model 
that are not discnssed in the paper. The waves are listed in Table D.l| together with their relative 
intensities as defined in Section IV C. Ont of the 69 waves, 64 have relative intensities below 1 %. 
The relative intensities of the 69 waves add np to a total of 26.4%. 


The employed analysis method has several known limitations. Some of them are related to 
presently open qnestions concerning the analysis of diffractively prodnced mnlti-body final states. 
Therefore, the systematic nncertainties indnced by the analysis method are in general difficnlt to 
qnantify. In Section IV F and Appendix B of the paper, we discnss systematic effects that affect 
the 18 selected partial waves. In the following, we focus on the discussion of systematic effects 
that are potentially affecting the remaining 69 waves. Since many of these waves have only small 
intensities, they are more susceptible to systematic effects. 


There are a number of systematic effects related to the truncation of the partial-wave expansion 
[see Eq. (24) in Section III C], i.e. to the choice of the wave set. There is currently no generally 
accepted objective method to determine from the data which partial waves actually contribute. 
As discussed in Section IV F, 17 of the 18 selected partial waves have been found to be fairly 
insensitive to changes of the wave set. This was assessed based on fits performed with a significantly 
smaller set of 53 waves. Mutual reshuffling of intensity (leakage) was observed, which, however, 
was shown to have no significant effect on 17 of the 18 selected waves. Typical examples of such 

and ID.6fl 


leakage effects are unphysical enhancements below 1 GeV/c^ as can be seen e.g. in Figs. D.3b 


Another class of artifacts that are related to the truncation of the partial-wave expansion are 
discontinuities due to the thresholds applied to 27 waves in the PWA model (see Table IX in 
Appendix A). These thresholds gradually reduce the waves set towards lower m^Ti for < 
1.7GeV/c^ and are required to stabilize the PWA fit (see also discussion below). 

Novel data-driven approaches to find an adequate wave set for a given data sample, like the one 
proposed by the authors of Ref. |66) . which is based on additional penalty terms in the likelihood 
function, might help to solve the issues discussed above and are currently under study. 

Non-resonant contributions, like e.g. the Deck effect, represent a continuum of waves including 
waves with very high spin. Thus, by truncating the partial-wave expansion, the non-resonant 
contributions are not fully accounted for. This might induce additional leakage. At the present 
state of the field, more theoretical work is required in order to treat the non-resonant contributions. 

The Bose symmetrization of the two final-state 7r“ breaks the orthogonality of the angular part 
of the decay amplitudes. In particular, amplitudes with the same quantum numbers but 

different decay channels may develop large overlap integrals, indicating that the corresponding 
phase-space distributions are similar. This typically leads to unstable maximum-likelihood 
estimates for the intensities of these waves. Such instabilities are most pronounced at low msTr, 
where only the low-mass tails of the isobars contribute to the decay amplitudes. Since the isobar 
amplitudes exhibit similar behavior at low it becomes difficult to distinguish them. In the 

current approach, this is counteracted by applying thresholds to some waves, thereby reducing 
the wave set in the region of low m^-K (see Table IX in Appendix A). 

The employed isobar model uses predefined parametrizations for the mass-dependent amplitudes 
of all isobars. Systematic effects due to possible imperfections in the isobar parametrizations are 
discussed in Sections IV F and VI and in Appendix B 3. Additional effects may be caused by 
small contributions from heavier isobars, like e.g. /9(1450), p(1700), or higher /q excitations, which 
are currently neglected. Including such excited states as isobars leads to unstable fits because 
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Table D.l: List of the 69 waves of the 88-wave PWA model that are not discussed in the paper. 
The intensities are evaluated as a sum over the 11 bins and are normalized to the total number of 
acceptance-corrected events. They do not include interference effects between the waves. 



Isobar 

L 

Relative 
intensity [%| 

Shown in 

0-+0+ 

p(770) 

P 

3.5 

Fig. 

D.la 


0-+0+ 

/2(1270) 

D 

0.2 

Fig. 

D.lb 


0-+0+ 

/o(1500) 

S 

0.1 

Fig. 

D.lc 


1++1+ 

[ttttJs 

P 

0.2 

Fig. 

D.2a 


1++0+ 

p(770) 

D 

0.9 

Fig. 

D.2b 


1++1+ 

p(770) 

D 

0.6 

Fig. 

Die 


1++1+ 

/o(980) 

P 

0.1 

Fig. 

D.3a 


1++1+ 

/2(1270) 

P 

0.5 

Fig. 

D.3b 


1++0+ 

/2(1270) 

F 

0.1 

Fig. 

Die 


1++0+ 

P3(1690) 

D 

0.1 

Fig. 

D.3d 


1++0+ 

P3(1690) 

G 

< 0.1 

Fig. 

d1^ 


1-+1+ 

p(770) 

P 

0.8 

Fig. 

D.4a 


2++2+ 

/2(1270) 

P 

< 0.1 

Fig. 

D.5a 


2++1+ 

P3(1690) 

D 

< 0.1 

Fig. 

D.5b 


2-+ 1+ 

Ms 

D 

0.4 

Fig. 

D.6a 


2-+0+ 

p(770) 

P 

3.8 

Fig. 

D.6b 


2-+ 1+ 

p(770) 

P 

3.3 

Fig. 

Die 


2-+ 2+ 

p(770) 

P 

0.2 

Fig. 

D.6d 


2-+ 1+ 

p(770) 

F 

0.3 

Fig. 

Dli" 


2-+ 2+ 

/2(1270) 

S 

0.1 

Fig. 

DlfT 

2-+ 1+ 

/2(1270) 

D 

0.2 

Fig. 

D.7a 


2-+ 2+ 

/2(1270) 

D 

0.1 

Fig. 

D.7b 


2-+0+ 

/2(1270) 

G 

0.1 

Fig. 

Die 


2-+0+ 

P3(1690) 

P 

0.2 

Fig. 

D.7d 


2-+ 1+ 

P3(1690) 

P 

0.1 

Fig. 

Dli' 


3++0+ 

M]s 

F 

0.2 

Fig. 

D.8a 


3++1+ 

M]s 

F 

0.3 

Fig. 

D.8b 


3++0+ 

p(770) 

D 

0.9 

Fig. 

Die 


3++1+ 

p(770) 

D 

1.0 

Fig. 

D.9a 


3++0+ 

p(770) 

G 

0.4 

Fig. 

D.9b 


3++1+ 

p(770) 

G 

0.1 

Fig. 

Die 


3++0+ 

/2(1270) 

P 

0.4 

Fig. 

D.9d 


3++1+ 

/2(1270) 

P 

0.4 

Fig. 

Dli" 


3++0+ 

P3(1690) 

S 

0.4 

Fig. 

Dlf] 


3++1+ 

P3(1690) 

S 

0.1 

Fig. 

D.9g 


3++0+ 

P3(1690) 

I 

< 0.1 

Fig. 

D.9h 



jpc jvre 

Isobar 

L 

Relative 
intensity [%] 

Shown in 

3“+ 1+ 

p(770) 

F 

0.1 

Fig. 

D.lOa 

3“+ 1+ 

/2(1270) 

D 

< 0.1 

Fig. 

D.lOb 

4++ 2+ 

p(770) 

G 

< 0.1 

Fig. 

D.lla 

4++ 2+ 

/2(1270) 

F 

< 0.1 

Fig. 

D.llb 

4++ 1+ 

P3(1690) 

D 

< 0.1 

Fig. 

D.lle 

4“+ 0+ 

Mis 

G 

0.3 

Fig. 

D.12a 

4“+ 0+ 

p(770) 

F 

1.0 

Fig. 

D.12b 

4“+ 1+ 

p(770) 

F 

0.4 

Fig. 

D.12e 

4“+ 0+ 

/2(1270) 

D 

0.3 

Fig. 

D.12d 

4“+ 1+ 

/2(1270) 

D 

0.1 

Fig. 

D.12e 

4“+ 0+ 

/2(1270) 

G 

< 0.1 

Fig. 

D.12ff 

5++ 0+ 

Mis 

H 

0.1 

Fig. 

D.13a 

5++ 1+ 

Mis 

H 

0.1 

Fig. 

D.13b 

5++ 0+ 

p(770) 

G 

0.3 

Fig. 

D.13e 

5++ 0+ 

/2(1270) 

F 

0.1 

Fig. 

D.14a 

5++ 1+ 

/2(1270) 

F 

0.1 

Fig. 

D.14b 

5++ 0+ 

/2(1270) 

H 

< 0.1 

Fig. 

D.14c 

5++ 0+ 

P3(1690) 

D 

< 0.1 

Fig. 

D.14d 

6++1+ 

p(770) 

I 

< 0.1 

Fig. 

D.15a 

6++1+ 

/2(1270) 

H 

< 0.1 

Fig. 

D.15b 

6-+0+ 

Mis 

I 

0.1 

Fig. 

D.16a 

6“+ 1+ 

Mis 

I 

0.1 

Fig. 

D.16b 

6-+0+ 

p(770) 

H 

0.7 

Fig. 

D.16c 

6“+ 1+ 

p(770) 

H 

0.1 

Fig. 

D.16d 

6-+0+ 

/2(1270) 

G 

0.1 

Fig. 

D.16e 

6-+0+ 

P3(1690) 

F 

< 0.1 

Fig. 

D.16f 

1++ 1- 

p(770) 

S 

0.3 

Fig. 

D.17a 

1-+0- 

p(770) 

P 

0.3 

Fig. 

D.17b 

1-+1- 

p(770) 

P 

0.7 

Fig. 

D.17e 

2++0- 

p(770) 

D 

0.3 

Fig. 

D.17d 

2++0- 

/2(1270) 

P 

0.2 

Fig. 

D.17e 

2++ 1- 

/2(1270) 

P 

0.3 

Fig. 

D.17ff 

2-+ 1- 

/2(1270) 

S 

0.2 

Fig. 

D.17g 


Intensity Sum 

26.4 
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the waves with the excited isobars have large overlap integrals with waves with the ground-state 
isobars, like e.g. /9(770) or [Trvr]^ (see also discussion of the /o(1500) isobar in Section IV B). 

Another effect that is neglected by the isobar model is a possible distortion of the mass-dependent 
isobar amplitudes due to final-state interactions. The above issues can be addressed in future 
analyses by extending the freed-isobar method that currently includes only isobars 

(see Section VI) to also include dominant partial waves with and isobars. 

First studies of the relativistic corrections to the decay amplitudes in the partial-wave analysis 
show that the 18 selected waves exhibit only small changes (see Sections III B and IV F). However, 
for some of the remaining 69 waves these corrections are not small and will have to be taken into 
account in future analyses. 


For completeness and for future reference, we show in the following for each of the 69 waves in 
Table D.l the intensity distribution summed over the 11 t' bins. The percent numbers given 
in the mass spectra are the relative intensities of the particular partial wave shown. In view of 
the systematic uncertainties possibly induced by the effects discussed above, we refrain from 
interpreting the intensity distributions here. 


D.l Waves with Positive Reflectivity 
D.1.1 = 0-+ Waves 



Figure D.l: The t'-summed intensities of partial waves with = 0 “*■ and positive reflectivity. 


D.1.2 = 1++ Waves 



Figure D.2: The t'-summed intensities of partial waves with = 1++ and positive reflectivity. 
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(a) 



Figure D.3: The t'-summed intensities of partial waves with 


and positive reflectivity. 
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Figure D.4: The t'-summed intensity of the partial wave with spin-exotic = 1 + and positive 
reflectivity. 
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D.1.4 = 2++ Waves 


2**2* n P 




Figure D.5: The ^'-summed intensities of partial waves with = 2++ and positive reflectivity. 


D.1.5 = 2-+ Waves 
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Figure D.6: The t'-summed intensities of partial waves with = 2 + and positive reflectivity. 
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Figure D.7: The t'-summed intensities of partial waves with = 2 and positive reflectivity. 
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Figure D.8: The t'-summed intensities of partial waves with = 3+“'' and positive reflectivity. 
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Figure D.9: The i'-summed intensities of partial waves with = S'*'’*' and positive reflectivity. 
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D.1.7 = 3 -+ Waves 
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Figure D.IO: The t'-summed intensities of partial waves with spin-exotic = 3 + and positive 
reflectivity. 


D.1.8 = 4++ Waves 



Figure D.ll: The t'-summed intensities of partial waves with = 4++ and positive reflectivity. 
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D.1.9 = 4 + Waves 



Figure D.12: The i'-summed intensities of partial waves with = 4 and positive reflectivity. 


D.1.10 Waves 



Figure D.13: The t'-summed intensities of partial waves with = 5'*"''" and positive reflectivity. 
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Figure D.14: The t'-summed intensities of partial waves with = 5+“'' and positive reflectivity. 


D.1.11 Waves 
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Figure D.15: The t'-summed intensities of partial waves with = 6+“'' and positive reflectivity. 












92 


D.1.12 = 6 + Waves 
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Figure D.16: The i'-summed intensities of partial waves with = 6 “*■ and positive reflectivity. 
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D.2 Waves with Negative Reflectivity 



2*^0-pinO) 7t D 2^+0"/^(1270)®P ^jg3 2++r/^(1270) ® P 



Figure D.17: The t'-summed intensities of partial waves with negative reflectivity. 
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